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In two earlier papers* in this series, the extent to which a square-integra- 
ble junction and its Fourier transform can be simultaneously concentrated 
in their respective domains was considered in detail. The present paper 
generalizes much of that work to functions of many variables. 

In treating the case of functions of two variables whose Fourier transforms 
vanish outside a circle in the two-dimensional frequency plane, we are led 
to consider the integral equation 



Ys? 



(x) = / J N (cxy) Vcxy<p(y)dy. 

Jo 



(*) 



It is shoivn that the solutions are also the bounded eigenfunctions of the 
differential equation 

a generalization of the equation for the prolate spheroidal wave functions. 
The functions <p {called "generalized prolate spheroidal functions") and the 
eigenvalues of both (i) and (ii) are studied in detail here, and both analytic 
and numerical results are presented. 

Other results include a general perturbation scheme for differential equa- 
tions and the reduction to two dimensions of the case of functions of D > 2 
variables restricted in frequency to the D sphere. 



* See Refs. 1 and 2. 

3009 
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I. INTRODUCTION 

Iii two earlier papers 1,2 in this series, the extent to which a square-in- 
tegrable function and its Fourier transform can be simultaneously con- 
centrated was considered in detail. In that analysis, the eigenfunctions 
and eigenvalues of the finite Fourier transform played a key role. These 
functions, denned for f re | ^ 1 by the integral equation 

MM = [/""hWv, (i) 

can be continued analytically throughout the complex plane. They 
possess a number of special properties that make them most useful for 
the study of bandlimited functions. The functions are complete in the 
class of bandlimited functions; they are orthogonal hi ( — 1,1) and also 
in ( — oo , oo ) ; the fy are also the eigenfunctions of the integral equation 

I. *<*-») (2) 

Zt 

whose kernel is positive definite; \f/ , the eigenf unction of (2) belonging 
to the largest eigenvalue, is in an appropriate sense most concentrated 
among bandlimited functions of given energy. These and other proper- 
ties are discussed in detail in Refs. 1 and 2. Some familiarity with these 
papers will be assumed in the following. 

In the present paper we consider certain aspects of the generalization 
of this earlier work to functions of many variables. Many of the structural 
results of Refs. 1 and 2 (as was pointed out there) depend only on the 
fact that the operator defined by the right of (2) is completely con- 
tinuous and positive definite. The generalizations to D dimensions are 
perfectly straightforward : we comment briefly on some of them in Sec- 
tion II, but do not belabor them. Our main concern here is with details 
of the explicit solution of some of the integral equations that are multi- 
dimensional generalizations of (1). An unexpected dividend of this work 
is that one of these equations is of interest in the theory of masers. 

In Section III, we point out some general features of the integral equa- 
tions to be considered. Section IV treats the case of functions of two 
variables whose Fourier transforms vanish outside a circle in the two- 
dimensional frequency plane. The analog of (1) is shown to be the in- 
tegral equation 
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w ( r ) = f J N {crr')y/crrW)dr\ 0<rgl. (3) 

Jo 

This integral equation also describes the modes in a maser interferometer 
with confocal spherical mirrors of circular cross section (Ref. 3, p. 488). 
The eigcnfunctions of (3) are shown to be the bounded solutions of 



. t ,, dV _ d*p . ( 22 i - N 2 \ _ 



(4) 



that vanish at x = 0. 

We call the solutions of (4) generalized prolate spheroidal functions. 
Section V is devoted to their study: 5.1 treats the case of small c; 5.2 
and 5.3 treat various asymptotic cases.* 

In Section VI, the results of Section V are used to discuss the eigen- 
values of (3). Various asymptotic forms for these quantities are derived. 

Section VII treats the case of functions of D > 2 variables whose 
Fourier transforms vanish outside a sphere in the D-dimensional fre- 
quency space. It is shown that this more general problem can be reduced 
completely to the case already treated in Sections IV, V and VI. 

Finally, in Section VIII we present some numerical detail about some 
of the eigenf unctions and eigenvalues encountered. Applications of these 
results will be presented elsewhere. ' 

II. GENERALIZATIONS OF EARLIER WORK 

We denote points in Euclidean space of D dimensions, E D , by vectors, 
x = {Xx,x 2 , •■■ ,.f D ). A square-integrable function of D variables, 
fix), is said to be K-limited if it can be represented as a Fourier integral 

fix) = (2tt) " / exp (ix-y)F(y)d!j (5) 

over the bounded region /?. Here xy = £ *#* is tne usual scalar 
product and we write dy for TJ r/ '/< • K / is of total energy A, then by 
Parseval's theorem 

A = f | fix) fdx = (2t)" d f | Fiy) \*dy, 

■'/■:„ ' • R 

whereas the energy off in the bounded region S is 

* Some of the results of Sections IV and V have been developed independently 
by J. C. Huertley," who was led to consider (3) from its laser applications. 
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f\m\ 2 dz 

= f dz(2w)~ 2D [ dxexp (iz-x)F(x) f dy exp (-iz-y)F(y) 

JS •> R Jr 

= (2tt)- c f dx f dy K a (x - y)F(x)F(y) 

J R J R 



where 



K, 



,(x - y) = (2tt)- d [ exp [iz- (x - y))dz 



(6) 



and an overbar denotes complex conjugate. The largest fraction of energy 
that an /^-limited function can have in the region S is therefore the 
maximum value of the fraction 

f dx f dy K a (x - y)F(x)F(y) / f | F(y) \ 2 dy 

J R •> R I •> R 

taken over all functions F square-integrable through R. This maximum 
is the largest eigenvalue of the integral equation 



W(x) = [ K s (x - y)Hy)dy, x € R, 

which is the analog of (2). 
The kernel (6) of (7) is positive definite, since 

I dx f dy K s (x - y)/(x)/(y) 

•> R J R 

= (2t)~ d I dz\ I dx exp (iz-x)f(x) 

JS \ J R 



(7) 



> 



whenever 



f l/(x) 

J R 



2 dx > 0. 



By well-known theorems (see Ref. 4, Chap. 6), the eigenvalues of (7) 
are real and positive and the eigenfunctions, orthogonal on R, are com- 
plete in the class of functions square-integrable in R. A complete dis- 
cussion of the simultaneous concentration of square-integrable functions 
in E n and their Fourier transforms can be given in terms of the largest 
eigenvalue of (7) as in Kef. 2, Theorem 2. 
The right member of (7) can be used to extend the domain of defini- 
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tion of \p. We define 

*(x) = i ( Ks(* - y)Hy)dy, x <E E D . 

A Jr 

Then for two different eigenfunctions of (7) 
[ *,(x)fc(x)d* 

Je d 

= — [ dx f dy ^(x)^-(y) / dz K fl (z - x)i? s (z - y). 

To evaluate the innermost integral here, we observe from (6) that K 8 is 
given as a Fourier transform, so that from Parseval's theorem, 

[ K s (z - x)K s (z - y)dz 

J E D 

= (2ir)~° [ du exp [-tu-(x - y)] = R s (x - y). 

Js 

One then finds 

[ J, t (x)Mx)dx = A- ( d* M*) [ dy A' s (x - y)fo(y) 

Jb d AiAj J K Jr 

-I [ dx^(x)fc(x). 

Ai J R 

The orthogonality of the fa over 7? thus implies orthogonality over E D 

as well. 

Other results of the one-dimensional case extend as easily to D di- 
mensions, but we do not dwell further here on this general structure. 

III. SYMMETRY CONSIDERATIONS 

In what follows, we shall be concerned with the explicit solution of a 
number of instances of (7). Considerable simplification occurs when the 
region R is symmetric, i.e., when x € R implies -x € R, and when S is a 
scaled version of R. We write S = cR where x £ cR if and only if x/c € R 
with c a positive constant. We restrict our attention henceforth to this 



case. 



Somewhat simpler than (7) is the integral equation 

a^(x) = / exp (icx-y)$(y)dy, x(fl (8) 

J R 
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which is a natural generalization of (1). We shall show in this section 
that solution of this equation is completely equivalent to solution of (7) 
when the symmetries just discussed maintain. We shall accordingly here- 
after take (8) as our equation of fundamental concern. 

From the symmetry of R, it readily follows that if ^(x) is a solution of 
(8), so also is \f/( — x), so that both ^,(x) = ^(x) + ^(— x) and 
^ (x) = ^(x) — \f/(— x) are solutions as well. The eigenfunctions of (8) 
can be chosen to be either even or odd functions of x. 

The complex conjugate of (8) is 



a\f/(x) = / exp (-icx-y)\f?(y)dy, x 6 R. 

J R 



(9) 



Multiply (8) by \f/(x) and integrate over R. Multiply (9) by ip(x) and 
hitegrate over R. Combining these equations, one finds on using the sym- 
metry of R that 

(a ± a) I \p(x)f(x)dx 
Jr 

= / dx / dy exp (t'cx-y)£(x)[<A(y) ± H~j)l 
J n j r 

If then \f/ is even, by choosing the negative sign in this equation, one 
obtains a — a = 0, whereas if \f/ is odd, by choosing the plus sign, one 
finds a + a = 0. The eigenvalues of (8) associated with even eigenfunctions 
are real: the eigenvalues of (8) associated with odd eigenfunctions are pure 
imaginary. It follows then that (8) is equivalent to the pair of equations 

Mc(x) = f cos cx-yf c (y)dy (10) 

Jr 

Mo{x) = \ sin cx-y\f/ (y)dy (11) 

Jr 

in which & and /3 are real. These equations have real symmetric kernels 
and we can fall back on the extensive theory in the literature treating 
such equations. We observe that the eigenfunctions of (10) must be 
even and that /3 C = cannot be an eigenvalue of this equation, for by 
Fourier theory the only even square-integrable function in R for which 

/ cos cx-y\f/(y)dy = 0, x £ R 

Jr 

is\f/(y) = 0. It follows then from the theorem on page 234 of Ref. 4 that 
the eigenfunctions of (10) are complete in the class of even functions 
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square-integral >le in R. A similar argument shows that the solutions of 
(11) are complete in the class of odd functions square-integrable in R. 
The solutions of (10) can be chosen real and orthogonal in R, as can the 
solutions of ( 11 ). Solutions of ( 10 ) are automatically orthogonal to solu- 
tions of (11) by symmetry. 

We have now shown that the solutions of (8) are complete in the class 
of functions square-integrable in R. The eigenfunctions can be chosen 
real, orthogonal, and either even ( in which case the eigenvalue a is real) 
or odd (in which case a is pure imaginary). We henceforth assume the 
\p so chosen. 

By iterating (8), one finds that the \p also satisfy 

X-Mx) = f /C(x - y)+(y)dy (12) 

•in 

with 

Kc ( x ) = QLY J exp (icz-x) dz = (2tt)- d J cr e itx dz (14) 

which is (7) in slightly altered notation and is the D-dimensional ana- 
log of (2). Since the solutions \f/ of (8) are complete, it follows that they 
arc also a complete set of solutions of (12). As was asserted, to solve (12), 
it suffices to solve (8). 

The eigenfunctions of (8) can be extended by demanding that equa- 
tion to hold for all x 6 E D . It is then easy to show that the extended \J/ 
are orthogonal in E D and that they are complete in the class of cfl-limited 
functions. 

IV. THE CASE D = 2, R A CIRCLE 

We now treat in detail the equation 

*Hxi,m) - ( e Mxm +* iV > ) +(yi,y2)dy 1 dy 2 (15) 

where R is the unit circle y? + y? ^ 1. Change to polar coordinates 
gives 



a+(r,e) - fdr'r' f dd' e ier " - ^'W/) 

.'o *'o 

= £r a w f d/ r'JAcrr') f* &f 9+" *(r'/) 



(16) 
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on making the usual Bessel function expansion. Here tfr(r,0) is exhibited 
as a Fourier series in 0. A simple argument then gives for the eigenfunc- 
tions of (15) and their corresponding eigenvalues 

&.»(r,0) = RoAr), <*o,n = 27rft,„ 

+N, n (r,e) = Rn,n(r) C ° S *j?2 <*„,„ = 2,rt*/V. (17) 

sin i v Uj 

N=l,2,..., n = 0,l,2,... 

where 



fanRxM = f JAcrr')R»Ar')r dr, = r = 1, 

Jo 



(18) 



n,N = 0,1,2,... 



All the eigenvalues of (15), except possibly the a , n have at least a two- 
fold degeneracy inherited from the symmetry of the circle. 
Our task now is to study the integral equation 

0R(r) = [ J N (crr')R(r')r dr , = r = 1. 

Jo 

It is convenient to make the substitutions 

7 - Vcfi, <p(r) = VrR(r) (19) 

to obtain the symmetric equation 

w (r) = I J N {crr') VcrVf (r) dr', = r = 1. (20) 

•>o 

Note that ip(0) = 0. We shall show that the eigenfunctions <PN,n(r) of 
(20) can be obtained as the solution of a Sturm-Liouville differential 
equation. 
Let 

K N {x) = J N (x)V^ (21) 

and let the operator M be defined by 

[Mf](x) = t K N (cxy)+(y)dy. 
Jo 

Denote by L x the differential operator 

d a 2^ d , (I ~ N* * A 

Lx = Tx {l ~ x) Tx + V~^~ ~ cx >- 
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Consider now 

+ fcj£ - cV)]*(y)dy 

= [K N (cxy)(\ - y-)+'(y) - cx(l - y 2 ) (22) 

■K„'(cxy)My)]U + /V(y) UV(1 - y 2 )K N "(cxy) 

- 2cxyK'(cxy) + f t ~ N - c 2 y 2 J KAcxy) dy 

where the right member is obtained by integration by parts. Here primes 
denote differentiation of the function in question with respect to its 
argument. The integrated expression vanishes if ^(0) = 0, since from 
(21), AVO) = 0. Also from (21) and the differential equation satisfied 
by Bessel functions, one has the identity 

K s "{cxy) = -(l + ^jp) KAcxy). (23) 

Substitute this expression in (22) to yield 

[MJ4\{x) = [ +(y)[-2cxyK'(cxy) 

+ (1 - N 2 + <W - cV - c 2 y 2 )K(cxy)]dy, (24) 
*(0) = 0. 
On the other hand, by direct calculation and use of (23), one has 

[UtyKx) = L x f K.Acxy)+(y)dy 

= I f(y) |"(1 - x 2 )c 2 ,fIW(cxy) - 2xcyK N '(cxy) 

+ (*-^ - o 2 J) KAcxy)~\ dy 
- f\(y) \-2*cyK K \*y) + |-(1 - x)c 2 y 2 
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= [ Hy)[-2cxyK N '(cxy) + (1 - N 2 + cxy 2 
Jo 

- ex - c\f)K N {cxy)\dy 

= [ML+](x) 

on comparison with (24). 

Let C be the class of functions square-integrable in (0,1) and twice 
differentiable there that vanish at the origin. Operating on functions in 
C, the operators M and L commute. It follows that solutions of 

L x <p(x) = —x<p(x) 

in C are also solutions of (20). Consequently, we next turn our attention 
to the differential equation. 



, 2s dtp d<p (\ — N 2 2 , \ n 



V. GENERALIZED PROLATE SPHEROIDAL FUNCTIONS 



(25) 



When N = \ \\\ (25), this equation reduces to the equation for pro- 
late spheroidal functions of order zero. We shall refer to bounded solu- 
tions of (25) for arbitrary values of N as generalized ■prolate spheroidal 
functions. These functions are similar in many respects to prolate sphe- 
roidal functions, as the development that follows shows. Bounded solu- 
tions of (25) exist only for discrete values of x, say xn.n , n = 0, 1, 2, . . . 
which we label so that xn.o ^ Xa\i ^ Xn,2 ^ .... We denote the corre- 
sponding eigenfunctions by ip N , n (x). 

5.1 Expansions in Powers of c 

Consider first the case when c = 0. Substitution of the series 



■ Say-" 



,a+2j 



into (25) shows that we must have a = \ ± N. If N ^ 0, the negative 
sign leads to solutions having a singularity at x = 0. If N = 0, a second 
solution can be found, but it has a logarithmic singularity at x = 0. 
We must have therefore 

a = \ + N. 
The coefficients are given by the recurrence 

(«+ 2j)(a+ 2;+ 1) - x 



dj+i = a 



1 (a + 2j + 2)(a + 2j + 1) + \ - TV 2 ' 



(27) 
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For large j, dj+i/aj —> 1, so unless the series terminates, this solution 
becomes unbounded as x — » 1 . Choosing x to terminate the series at 
x a+2n , we have 

X = X.v,„(0) = (N + 2n + J) (N + 2n + |) (26) 

for the eigenvalues of (25) when c = 0. The series solution now becomes 
(when a„ is set equal to unity)* 

V = 7\., n (.T) = .r v+! 7?.v. n (.T) 

HsA-v) = F(-n,n+N + 1;JV + 1; x 2 ) 

where 

PM v , , afc z , fl(a+ Db(b+ 1) 2 2 , 

is the usual Gaussian hypergeometric function. The polynomial R N , n (x) is 
readily expressed in terms of Jacobi polynomials P„ c (a:) (Ref. 5, 
Chap. IV). Adopting the notation of Szego, we have 

R*M) =( n ^) _1 Pn W ' 0) (l -2.r 2 ). (28) 

From (27), (28) and the known properties of the Jacobi polynomials, 
one finds 

7V,(1)= (-!)"(* t^r (29) 



J » Of- 



8bb' 



2(2w + JV + 1) 



(■ : N ) 



(30) 



2(n + AT+ l) 2 (2n + A r )/?.v. n+1 

= (2n + -V + 1 ) [(2n + A< ) (2n + AT + 2) (1 - 2.r 2 ) 
+ A^]i?.v.„- 2n 2 (2« + N + 2)/r\ v ,„_ 1 

(2n + A^).i:(l - x 2 ) ^-R.s:„(x) 
ax 

= n[(2n + A0(1 - 2a; 2 ) - iV]ft WfB (x) - 2n 2 R v .„-,(.r) 

.r 2 r.v.„(a-) = 7-v, n I 7 , .v. n+ ,(.r) + 7 .v,>.v.„(.r) + y^Ts.n-M (31) 

* It has been called to our attention that our T N , n (x) are closely related to the 
Zernike polynomials. These latter arise in the diffraction theory of aberrations. 12 
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(n + JV + 1) 

lN,n = — 



(2n + N + l)(2n + JV + 2) 
y*.„° = 2 V 1 + (2n + AT)(2n + AT + 2)) (32) 



-, 

7AT,n = 



(2w + AT)(2n + N + 1) 

I T N , n (x) | £ 1 for ^ x ^ 1. 

The function T Ntn (x) has n zeros in (0,1]. We define T N , n {x) = if 
n < 0. 

Returning now to (25) for arbitrary values of c, we attempt a power 
series solution in c by writing 

<p(x) = *,„(*) - ZV.„(aO + Ilc 2i Qj(N,n,x) (33) 

00 

X = X*.»(c) = x*..(0) + Ec'W^.n), (34) 

where the Q's and a's are independent of c. When this latter quantity is 
zero, this solution reduces to that already found. As is shown in Appendix 
A, the Q's and a's can be determined recursively in an elementary man- 
ner. We have 

Qj(N,n,x) = £ V(tf,tt)7V»+*(a:) (35) 

with 

aj(N,n) = £ AV'-'tAr.nh^n-/', .7 = 1,2,... (36) 

\XH,n+m(0) - x»A0)]Aj(N,n) 

= iTa k (N,n)A m ^ k (N,n) - £ ^_ fr+m ^ I (Ar,n) 7w .„-. +m fc , (37) 

fc-l A=-l 

w = -ii -i + 1, . . . ,i; j = i, 2, . . . . 

Here A k 3 (N,n) is defined to be zero if | fc | > j, or k < — n or k = 
and j ?* 0. In addition we have A m °(N,n) = 0, m ?± 0, A °(N,n) = 1, 
a (AT,n) = 0. For use in (36) and (37), the 7's of (32) must be defined 
so that for n < 0, y N , n l = yn.n = 7jv.n+f * = 0. 

To terms of order c the eigenf unctions and eigenvalues of (25) are 
explicitly 
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XA ,„( C ) - Un + N + 1) (2n + N + |) 

2 (38) 

+ K 1+ (2n+A0cl + iV + 2) ) c2 + Q(c4) 

/ v _ r / v , / ?t 2 r A -,„-i(ar) 

W - W ~ 1n - {X) + V4(2n + iN0 f (2n+JV + 1) 

(n + AT + l) 2 2y„ +1 (a0 \ c 2 + ^ 



4(2n + N + 1) (2/i + iV + 2) 2 
In view of (35), the series (33) can be formally regrouped to give 

«r,(«) = Ed/' n (c)^.,(x). (40) 

o 

Substitution in (25) yields the three-term recurrence 

2 1 j iV.n 

CjN.j-l Uj-l 

+ [A*,/ + (2j + AT + 1) (2j + JV + |) - X ] ch N ' n (41) 

■+- cyjfj+i dj+i = 0. 

This recurrence can be used to determine the dj N,n and the eigenvalues 
in a manner quite parallel to that used in the study of prolate spheroidal 
wave functions. The method of Bouwkamp 6 can be adopted and used 
advantageously for the computation of the dj '" and the eigenvalues for 
values of c too large to permit effective use of (33) and (34). The r/'s 
can, of course, be expressed in terms of the A'b of (35). One has 

d«+/ ,n (c) = f) A/(N,n)c 21 , j = -n,-n+l,.... (42) 

The series solutions (40) or (33), (35) for the generalized prolate 
spheroidal function are, of course, valid only for ^ x ^ 1. To obtain 
a series valid for x > 1, we use (20) and the fact (established in Appendix 
B) that 

/' JAc,y) V^y TV., («)* = ( N + "V ^±^ . (43) 

The solution (40) then extends for all x by the series 



IpH.n U) = 2_, Oj 



C70* 



(44) 
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which is obtained by inserting (40) in the right of (20) and integrating 
term by term. 

The eigenvalue y N ,n can be expressed in terms of the dj '". Divide 
both sides of the equation 

yN.n(pft,n(x) = I Jn(cxx) y/cxx' (p N , n (x)dx (45) 



by x N+i and take limits as x — > 0. From (27) and (40), we see that the 
left member of (45) becomes 

7».» l>/'"fi".n(0) = 7".* E <*/'"• 
Since J«(x)Vx ~ (x/2)"\/x/r(N + 1), the right of (45) becomes 



JV+* 

C 



T(N 4- 1)2* p 



-£e*/' n f T N , (x')T N j{x')dx 

" ,=0 ^0 



r(JV + 2)2" +1 

where we have used successively (40), (27) and (30). The combined 
result is 

c N+i d N - n 

y». n = 2 " +1 r(Ar + 2) £ <*/•" ' (46) 

The integral equation (45) is also useful for obtaining the asymptotic 
behavior of <p N , n {x) as x — * °o. We have 

y*a0**\V = — duu j nW uzi — (47) 

on letting era;' = w. Now (w w+1 Jy+i)' = (u N+1 J N ), so that (47) can be 
integrated by parts to yield 






<PN,n (u/cx) 



iff+i 



'Jo duu Jn+m Tu sea J 
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= — y/cx <pN, n (l)J n+\(cx) - R. 
ex 

For large x, this becomes 

, v ,-v , /2 cos [cz - (N + l)(r/2) - (tt/4)] 

yit,n<pN,n\.X) = tpN,n\l) A/ ~ 

(48) 



»&• 



This of course is consistent with (44). If now we define <p N , n (x) to be a 
generalized prolate spheroidal function normalized so that for large x 

„ . cos [ex - {N + 1)(t/2) - (>/4)] ,. q * 



(48) gives us 



c.r 



Tv V?^.n*(D (50) 



a relation that will be useful to us later. 

5.2 Asymptotics for Fixed n and Large c 

The behavior of generalized prolate spheroidal functions for large c 
can be determined by methods quite parallel to those used in Ref. 7 in 
discussing the prolate spheroidal functions. Five different asymptotic 
forms for <p N , n {x) are found, depending on the x range under considera- 
tion. These are properly joined to furnish a solution for all x. For most of 
these regions, we content ourselves here with writing only the leading 
term of the asymptotic development. 

In (25) we make the substitution t = xy/c. There results 

U - (l/c)M* + (x/c)<p = (51) 



where the operators are given by 






L = 4r, + *-r~ - ? 



Now the equation 



*-'£+*£■ 



LU + \U = 



(52) 
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has solutions 

U = UNin (t) = e^'Y^L^it 2 ) 

(53) 
X„ = 4n + 2N 4- 2, n = 0,1,2, ... 

(see Ref. 5, p. 99) where L n (x) is the Laguerre polynomial of degree n 
in Szego's notation. The function Us,n{t) has n zeros in (0,°o). This 
suggests attempting solution of (51 ) for large c by the series 

<p N ,n = U N , n (t) + £ (l/c y )£y(iV>,0 (54) 

3=1 

Xw .n(c)/c = 4n + 2JV + 2 + £ (l/c%(A/». (55) 

We now note that 

MU N ,n(t) = MA\« Us.n+1 + MAT.n £^V,n + AtJV.n" U N ,„- 2 (56) 

where 

M.V,,, 1 = (w + 1) {n + 2) 

Wr.. - = -[(2n + 1) (n 4- N + |) + |] (57) 

ns,n~ l = (n + N) (n + N - 1), 

a fact which can be readily derived from (52), (53) and the properties 
of Laguerre polynomials. The perturbation scheme of Appendix A applies 
therefore, and we find at once that 

Sj(N,n,t) = T,B k j (N,n)U N , n+2k (t) (58) 

k — j 

where the J5's and fr's are given by the recurrence 

bi(N,n) = £ B-S~ l (Nfi)pm0-» k , j = 1,2, ... 



SmBj(N,n) = £ b k (N,n)Bj- k (N,n) 

- 2 B_ k+ „r l (N,n)» N , n+2 k (m - k) (59) 



i 

r 
™ = -ii — i + i i; i = i, 2, . . . 



with the convention B k \N,n) s if | fc | > j, or fc<— norfc = 
andj * 0. We take B m °(N,n) = 0, m ^ 0, B (N,n) = 1, b (N,n) = 0. 
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In this manner we obtain explicitly 
Xf ,, n (c) = (4w + 2N + 2)c - [(2n + l)(n + N + I) + f] 

_ (A/ + 2w + l)[2n 2 + 2rc(AT + 1) + N + 2] + Q A \ (60) 

which gives the behavior of x.v, n for large 6-. 
We write the solution just found as 

;=1 C t=-j (hi) 

The right side of (61 ) is ordered in powers of c _1 when expressed in terms 
of the variable /. However, if / = x\fc is substituted, the terms are no 
longer so ordered since Ut/,m+i(x-y/c)/U tr ,m(z'\/c) = 0(c). The range 
of x values for which the first few terms of (61) furnish information 
about ip N „ vanishes as c gets large. We shall use (61) only for ^ x ^ 
l/c>. 

To obtain an asymptotic form for <p Ni „(x) for c ^ x ^ 1 — (1/c) 
it is convenient to write <p N , n (x) = x N+1 \f/ N§n (x) and set y = y/\ — x 2 . 
Equation (25) now becomes 



(i-/)g + [l-(«r + .),]g 

+ \x».n - C - (N + |)(iV + |) + 6-V] * - 0. 



Into this equation, substitute x*.»(c) as given by (60) and set 



g*(l - y)" 



One finds then for t>, 



Accordingly we write 



dy 



+ o(D=o. 



,,-v+S(l - /y )»^ 

**•■ (x) ~ V»(i + z/) ,v+n+1 
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To obtain an asymptotic form for <p N ,„(x) valid near x = 1, set 
y = s/c in (62) and again use (60) for x*,n . There results 

Accordingly we write 

tp N ,n\x) ~ x N+h I (cy) 

d (64) 

y = Vl - x\ 1 - (1/c) ^ a; ^ 1 

where /„(x) is the modified Bessel function. (See Ref. 8, Vol. II, p. 5). 

When x > 1 , we set z = \/x 2 — 1, and have ?/ = iz. The solutions 

<Pa',„ 2 and ipN.n then give rise to two more asymptotic forms. We write 

<PN.n(x) = x N+i J (cz), 1 ^ x £ 1 + - c (65) 

"-■<«> = ^ Re v^iV^ : - +1 ■ ' + 1 * x > (66) 



2 



= -v/aJ^TL. 



We now determine the joining factors for these five solutions. In 
<PN,n and ifN.n we set x = u/& and let c become large for fixed w. One 
finds 

1/ / l\ \ — l-)C 2n+N+l -i/2y^/2 

^#,, U/c 1 ) ~ j u *e 

c-(2n+N+&H 

^, n (.c = M /c ) 2N+2n+l u e 

where we have used the fact that 

L n s {uW~c) ~ (-llVV/nl 
When /y = v/\/c, one finds for fixed v and large c 
<pN,n(v = v/y/c) ~ cV B /V» 

¥>*.,< "'(z/ = y/Vc) ~ /^-i e* /\/v 

where to obtain this last expression we have used the known asymptotic 
formula I„(x) ~ e x /\/2^x (see Ref. 8, Vol. II, p. 86). Finally, when 
z = v/\/c we find 
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ips, n \z = v/y/c) ~ \/ - 



2 -4 cos (vy/c — 7r/4) 
— c — 



<Pn 



,„ a (z = v/y/c) ~ c 



y/v 
l cos (v y/c — 7r/4) 
y/v 



where we have made use of the formula (see Ref. 8, Vol. II, p. 85 i 
J (z) ~ («/2) H cos (2 - 7T,. 4). 

All these results can he summarized in the following statement: 



0N,„(X) 



-|S/2.,V + J r (N)/,2\ 

e t L„ (.t ), 
x N+ k c "(l - >/)" 

h Vj(i + v) N+n+] 

k 3 x N+i I (cy), 

k<X N *J (cz), 



^ x ^ c ' 



c ^ a ^ 1 — c~ 



1 - c" 1 ^ a; ^ 1 



/.V v,} He -7= 



e'«(l - &)■ 



i*(l + IZ) 



I +c~ l £ X 



1 ^ x ^ 1 + c 1 
(67) 



where 



/ = xy/c, y = y/\ — x 2 , 

( _ 1) „ 2 . v+2 , (+lcn+ .v/ 2+ i e -, 



= y/x 2 - 1 



/-••> = 



til 



, , (-i)V^ +2n+3/ v + " /2+i <r 

/,;, - fc< - - - 



u = 



(-D"2- v+2n+ y +im+ v 



n : 



is the asymptotic form for large c of a bounded continuous solution of 
(25) belonging to the eigenvalue (60). 
We next calculate the normalization constant 



1 



N N " 



[0lf,n(x)f dx. 



For the contribution due to ^>. Vi „ we find 
_i ■ 

h V c Jo 

= VT rdtic-"u K [L n (N) (n)) 2 
2y/c Jo 

= r(» + N + i) (1 + 0(c » +v v<) 

2y/cT(n + 1) 
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where we have used the fact that 

r x „,, «/ u z , r(n + « + 1) 

/ c x [Ln {x)\ ax = — — — -r — - 

Jo T(n + I) 

(see Ref. 5, p. 99). It is not hard to show that the contribution to 1/Ns,n 
from integration over the region c 4 ^ x ^ 1 is 0(c p e _v ^) for some p > 
0. We have then 

,_ 2Vcr(n+ 1) rfiR * 

5.3 Asymptotics for n and c Both Large 

The techniques employed here again follow very closely those used 
in Ref. 7. We accordingly give a minimum of detail. 
We assume that when n and c are both large x can be written 

xs.n ~ c + 25c + b + fo/c + • • • . (69) 

The ranges of n and c for which this is valid will appear in the analysis 
to follow. 

In (25) make the substitution x = t/c and replace x by (69). One 
finds 

and hence for large c, <p(t) ~ -\/~tJrf(t). We write 

¥>*.» 6 (.r) = VxJAcx), 0£x£-7=. (70) 

Returning to (25) with x replaced by (69), we observe that the sub- 
stitution 



«p[<(--|^f+i)] 

f / 



v~ 



yields — -f ( - I = 0, so that for large c, y becomes constant. After 
ax \c / 

multiplying this solution by a complex constant, we take its real part 

for the next section of <p. Explicitly we define 

,, , . / 2 cos [ cr ~ I log rri ~ (N + * } l] (7D 
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Note that when x = u/y/c and c is large (71) becomes 

The asymptotic formula for J N (see Ref. 8, Vol. II, p. 85) shows that 
*»..'(* = u/yfc) ~ ]/^ c cos [ w Vc - (N + J) |] 

also, so that v'y... 6 and tp s ,n agree for large c in the neighborhood of 

x = l/y/'c. 

To find an appropriate asymptotic form for <p valid near x = 1, 
substitute tp = .v s ' +, c icn ~ x) u into (25) with x given by (69). Now make 
the substitution x = 1 — i£/2e. There results 

Accordingly, we are led to define 

.»(*) = *» + V' M *[j - *|, 1; -2fc(l - *)] 



where 



(72) 



, . oi , o(o+l) 3 2 , 

<Ha,lr y r) - 1 + £ fl + &(& + 1) 2! + 



is the confluent hypergeometric function in the notation of Ref. 8, Vol. 

I, Chap. 6. 

The solution (72) is real. Its asymptotic form for large c when x = 1 ± 
v/y/c can be found from the known f asymptotics for the <?> function. 
One finds 

/H ±5(t/4) 

8/ , i /~\ V2e 

^ ( * = 1±v/ ^ c) ~V^W) (73) 

cos \vyfc =T= | log (2»V3) ± «(«) - 5 J 
where the real functions R(S) and 0(5) are defined by 

r Q + i-0 = fi(5)6* tf \ (74) 

This latter definition is made precise by requiring 0(5) to be continuous 
with 0(0) = 0. 

Now when x = 1 - y/Vc. I" 1 ) shows that 

t See Ref. 8, Vol. I, p. 278, Eq. (2). 
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<p N , n \x = v/Vc) ~ J^ 



COS 



vy/c + ^ log - 



(N 



»H- 



(75) 



2 *"° 2y/~c 

Comparison of this expression with (73) shows that <p N ,n and ( — 1)*- 
R(8)e s * ^(pN , n s /2\/ ir are asymptotically the same for x = I — v/y/c 
provided 

c + 5 log(2V^) - 8(8) = (N + 1)(tt/2) + *« (76) 

with 5 an integer to be determined shortly. 
Quite analogous to (71) is the solution for x > 1, 

e ,(W» cos [ex - I log ^i - (JNT + 1) I - l] 

Vrc y/x 2 - 1 

When x = 1 + v/y/c and c is large, this solution becomes 

g5(ir/2) 



(77) 



oos[c+»V3-|log^-(Ar+l)|-|]. 



Comparison with (73) shows that this is the same as ( - 1 )"(R(8)e s " li /2 • 
y/Tr)<PN, n % (x = 1 -f- v/y/c) when account is taken of (7fi). 
Our results thus far can be summarized as follows: 

(y/xJAcx), ^ x ^ c~ h 



-cos 



r 5 , I - x 



- CV + » 5 



I] 



7TC 



vr^" 



C * ^ .T < 1 — C 



,-I 



(-l)"R(d)e 



,S(x/4) 



V.V 



',•(*) ~S 



2Vtt 



J/H-.gfctt-.) 



*[l-t| f l;-2fc(l-*)] f 



(78) 



.T - 1 



,.-; 



""» cos [ - - 1 log r+i - (Ar + :) I _ f] 



V^rc 



X + 1 



a; ^ 1 + c 



-i 
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is the asymptotic form for large n and c of a continuous solution of 
(25) provided 5 and q are chosen to satisfy (76) and the requirement 
that ip as given by (78) has n zeros in the open .r-interval (0,1). The 
corresponding eigenvalue is given by x.v.n ~ c" + 25c + 0(1). Higher- 
order terms can be found by methods analogous to those presented in 
Ref. 7. 

When c becomes large and 5 remains fixed, i.e., 5 = 0(1), the number 
of zeros of <p N , n ( x ) m ° < x = 1 can be estimated roughly from (78). 
Using the asymptotic expansion for J N , we find that <^,„ 6 (.t) contributes 

z* = (V~c/*) +0(1) 

zeros as x ranges from zero to \/-\/c. From <pjv.„ (x) we find 

Zl = (1/ T )[ C - 2V~c + (5/2) log Vc] + 0(1) 

zeros for \/y/c ^ x ^ 1 - 1/Vc- Finally, by using the asymptotic 
form (73) for <pat,„ 8 , the number of zeros of <p for 1 — l/\/c ^ x ^ 1 
is estimated as 

z 8 = (lM[V~c + (5/2) log Vc] + 0(1). 

Since we must have n = z 6 + z 7 + z s , the last three equations show 
that 

nir = c + 5 log 2\^ + 0(1). 

Combined with (76) this implies that as c — > <», 

0(5) + (iV + 1)(tt/2) + irq - n-K = 0(1). (79) 

The eciuation just established can be used to obtain a limiting result. 
Let N be fixed and suppose that n grows with c according to 

n = (l/ir)[c + frlog (2Vc)] (80) 

where b is a fixed number (independent of c). Multiply this equation 
by 7r, add to (76) and rearrange to obtain 

(5 - b) log (2-\/c) = 0(S) + (#+ 1)(t/2) + tt 9 - mtt = 0(1) (81) 
where the last equality comes from (79). Divide (81) by log(2-\/c). We 
then obtain the limit result: if n grows with c according to (80), then 

lim 5 = b. (82) 

C-*30 

VI. ASYMPTOTICS OF y N ,„ AND \ N ,„ 

6.1 Fixed N and n. Large c 

The asymptotic solution <p Nt Jx) given in (67) has the values 
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0N,n(i) = />'3 

X 

On recalling the definition given in (49), we see that for large c 

<Pn,n*{x) ~ (-l)"^iv.„(a:)/cfcB, 
so that for fixed n and iV as c becomes large 

^•U)~^-C-l>-4/£. (83) 



Equation (50) then gives 



( — 1Y 

y N , n ~ ^^ . (84) 

Vc 



We now proceed to use (84) and the useful formula (to be established) 

^«^[*r, f (l) " 1] (8.5) 

ac 2c 



where 



f <p N , n \x) dx = 1 (86) 

Jo 



to get a much stronger statement regarding the asymptotic behavior 
of y N , n . First we establish (85)-(86). 

For simplicity of notation let us write (45) as 

7n<P,Xx) = [ K(cxx')<p„(x') dx' (87) 

Jo 

where we have suppressed dependences on N. Differentiating, we find 

— <Pn{X) + 7„ 

dc ac 

i , (88) 

- t xx'K'{cxx')<p n {x')dx' + [ K(cxx')^4^ dx '- 
Jo Jo ac 

Differentiating (87) with respect to x gives 

y n <Pn(x) =- [ xx'K'(cxx')<p n (x')dx', 

x Jo 

so that (88) becomes 
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dy„ / A , d<P,X-r) -r ,, v , f »- /„__/% d<p„(x') ,. 

— - <p„{x) + y„ — - — = - y,,<p„ (x) + / K {cxx ) — dx . 

dc oc c Jo oc 

Multiply this equation by <p„{x) and integrate. One finds 
oc Jo Jo oc 

.'o 2c o.r Jo oc 

where the last term has been obtained by interchange of orders of inte- 
gration and use of (87). Equation (85) then follows by integrating the 
first term on the right by parts and by using (86). 

To use effectively (85)-(86) it is convenient to introduce kn,„ = ( — 1)"- 
\/cyN, n . We then have 

J_**HL= * ,,„»(!) (89) 

K N , n dC 2C 

lim k n .„ = 1 (90) 

from (85) and (84) respectively. From (67) and (68) we see that 

n 2A4-4n+4 AT+2n+2 -2c 

2/1 > 7. 2 AT 2 - T2 C e 

<P„,n(D~foN„, n - r(B+1)r(n + iV+1) - 

Using this expression in (89) and integrating, we obtain 

. |00 7rJ / .N+2n+\-2t j. 

log ""■" '• = r(»+i)r(» + y + i) /. ' e * 

Integrating by parts and using (90), we finally find 
(_!)» (_i)» 7r 2 2W+4B+2 C Ar+2n+i e - 2c 



5H© 



(91) 



7 "" v^ r(n + l)r(n+tf + 

In terms of X of (13), we find from (17) and (19) 

X*.„ = cys.n (92) 

so that 

■,2A'+4n+3 W+2n+l -2c 



\ - 1 - T 

A/V.n — 1 



r(n + l)r(n + N 
6.2 Fired JV and n. £maM c. 



^['«©: 



(93) 



We use (46) to obtain an expression for y Ni „ for small c. From (42), 
it follows that 
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d N ' n (c) = Z Aj(N,n)c 2! = A- n n (N,n)c 2n [l + 0(c 2 )] 

(94) 
^ (-l)T(n + l)T(n + N + l)r(JV + 2) 2n[ , ,., 

2 2 T(2n + AT+ l)r(2n + iV + 2) L "^ V n 

where we have used (124), (26) and (32). From (42) one has 

Z d/'" = d/'" + 0(c 2 ) = 1 + 0(c 2 ). (95) 

Equations (94), (95) and (46) now give 

(-i)"r(n + i)r(n + W + i) 2 „ + ,v-H 



TAT.r. = 



2 2 " +JV+1 r(2n + AT 4- l)r(2n + N + 2) 



c 2 " +K+ \ (96) 



Higher-order terms could be obtained in a similar manner. An alterna- 
tive route, however, is to use (85) and (86). From (39) and (30), one 

sees that [2(2n 4- AT + l)] 4 ( n + N J (- 1)"[1 + 0(c*)] is the normaliza- 
tion factor for (39). Using (39) one then finds for a normalized solution 
PAr.nU) = (-l)V2(2n + N + 1) 

L 1 + 4(2w + N)\2n + AT + 2) 2 J + ° (c )- 
Inserting this expression in (85) and integrating, we find 

= (-l)T(n + l)r(n + N + l)c 2n+iy+ ' 

7*, n 2 2n +' v + 1 r(2n 4- N + l)r(2n 4- A^ + 2) 

(97) 

. 4(2n + Nf{2n 4- tf 4- 2) 2 "•" V '. 

6.3 Asymplotics for n and c Both Large 

To obtain an expression for y N ,„ valid for n and c both large, we use 
(77) and (49-50). For the asymptotic solution (77) we have 

and for very large x 

e Hrm cos [ cx _ (j y + l)( x /2) - (tt/4)] 
Vttc .r 

Comparison with (49) shows that <p N ,n = s/ir/ce~ w tp NlJt , and (98) 
and (50) now give 
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(-!)«#( 8) e-« ( " 4 > 



iK.n 



V2tc 



(99) 



Now (91) and (97) show that for large and small c the sign of y N , n 
is the same as the sign of ( — 1)". As c varies, 7a-,„ cannot change sign, 
for by (92) if y N ,„ were to vanish for some value of c ^ 0, so would 
X A -,„ . Since, as we have noted in Sections II and III, the kernel K c of 
(12) is positive definite, this is impossible. We can therefore replace 
q by n in (99) and we have 

q = ;t(mod2). (100) 

From the definition (74) of R (5), one has 

[*Wr=r(| + i0r(I- i i) = S5r ^. (ion 

Here we have used the functional relation [Ref. 8, Vol. I, p. 3, Eq. (7)] 
for the gamma function 

y(\ + z)Y{h - z) = ir sec 71-2. 

Equations (99), (100), (101) and (92) combined are 

Tat,,, ~ / t= , =f f , X/v.n ~ i j ?i. (102) 

Finally from (80), (82) and (102) we have the limiting result: if 



n 



= [(l/r)(c + blog2Vc)] 



where the brackets denote "largest integer in" and b is a fixed number, 
then 

^"•""rqh 5 " (103) 

VII. THE CASE D > 2, R THE UNIT SPHERE 

In the previous sections, we have treated the important special case 
D = 2, R the unit circle in considerable detail. Most of the analysis there 
was concerned with solving the integral equation (20). Fortunately, as 
we shall now see, the solution of that equation also affords a complete 
solution of the case R the unit sphere centered at the origin in E D , 
D — 3, 4, .... In treating this general case, we shall draw freely on the 
theory of D-dimensional spherical harmonics as given, for example, in 
Ref. 8, Vol. II, Chap. XI. We follow the notation of this work and set 

D = p + 2, p = 1,2, ... . (104) 
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Let x = r£ and y = r'n where I and n are unit vectors in E p+2 . Equa- 
tion (8) now becomes 

af(rjO = [ dr'r ,p+1 [ exp (icrr'Z-n)+(r', n ) dtt(n) (105) 

where ft is the surface of the unit sphere in E p+2 . 
Now let 

h(N,p) = (2V + p) (N + f N ~ 1)! , AT = 0,1,2,..., (106) 

and let Sh (£)> I = 1,2, ■ • • , h(N,p), be a complete set of orthonormal 
surface harmonics of degree N. The Funk-Hecke theorem (Ref. 8, Vol. 
II, pp. 247-248) asserts that 

f exp (icrr'Z -n)S N l {n) dSl(n) = H N (crr')S N l (t) (107) 

Jn 



where 
H N (crr') 

2 7r (p+D/2A/!( p _ i)i 



m 



(JV + p - 1)! 



f e iCTr ' u C N pl \u){\ - u 2 ) {p - im du (108) 



is independent of I and CV(w) is a Gegenbauer polynomial (Ref. 8, Vol. 
II, p. 235). By expanding \p in surface harmonics, 

at=o «=o 

we find from (105) and (107) 

«».i«#.i(r) = f drV^^CcrrO^y.^r') , (109) 

■>o 

from which it is seen that Rn.i(t) and a*,* are independent of /. We have 
the expected degeneracy of eigenvalues due to spherical symmetry. 
Now [Ref. 8, Vol. II, p. 236, Eq. (25)] 

tN {u) - ~^^7p + l\ ¥ , U } d^ K ; 

where (a) N = a(a +1) ••• (a -\- N — l),so that from (108) 

(p+l)/2/ 1 \N ,-1 iN 

ZlT \ — I) / j icrr'u Ci /, . ,2\ A r +(p-l)/2 



ffAr(crr') = ,.. , ,, v , \ 1X / d« e ,err u |^ (l - « z ) i 



r(E±l)^(^l)^ 
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Integration by parts gives for the integral here 

(-1)* C du(l -u l ) NMp - l)l2 ~e icrr ' u 
J-i dn N 

= (-icrrT f due icrr ' u {\ - uY +(v - im 

= (-i)Vrf (n + P±i) 2 w+p/2 (crr')- p/ V^ p/2 ( C rr') 

where we have used the Poisson formula 

r(f + \)J v {z) = ir"*(2/2)' /_ e" u (l - a 2 )*" 1 rfu 

[Ref. 8, Vol. II, p. 81, Eq. (7)]. We have then, finally 

H N (crr') = i s \2ir) x+ "' 1 Js +v „{crr')/{err') pi \ 

We see now from (10 ( J) that the eigenfunctions and eigenvalues of 
(105) are 

4>*,i,n(j£) = Biv.»(r)S JV '(?), / = 1, 2, ... , h(N,p) 

a N ,„ = i N (2*) 1+pl %. n (110) 

N,n = 0, 1,2, ... 
where 

fl*.JtM.M = [ l j ^ P ' 2 ^ ) r' p+l R N , n (r') dr'. (Ill) 

Jo (err) pl - 

These equations are the analogues of (16), (17) and (18). Set 

7 = /3c (p+n/2 , * = r lv+1), *R. (112) 

Equation (111) becomes 

W (r) = f .7 .v + ^(ar')\/c7^(f') dr'. (113) 

This, however, is (20) with .V replaced by N + p/2. The formulae of 
Section IV for the solutions of (20) can be taken over exactly replacing 
AT by N + p/2 throughout. (Expressions involving factorials must be 
replaced by the appropriate ones in terms of T functions when p is an 
odd integer.) Together with (110), (111) and (112), they provide solu- 
tion of (10.5) for all /) ^ 2. 

It is interesting to note that the one-dimensional case treated in Refs. 
1 and 2 can be obtained as a special case of the present theory by ap- 
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propriate interpretation. The parameter N of this section is the degree 
of the homogeneous polynomial solution to Laplace's equation in D 
dimension afforded by the spherical harmonic Sn when expressed in 
rectangular coordinates. When D = 1, Laplace's equation d\J//dx = 
has only two homogeneous solutions, \f/ = k and \J/ = x, respectively of 
degrees zero and one. For D = 1, i.e., p = — 1 from (104), we have 
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Fig. 1 — Curves of xtf.n of (25) vs c. 
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only two allowed values, N = and N = 1. The quantity N + p/2 
occurring in (113) then has values -\ and \. The kernel becomes 
y/2jr cos err' and -v/2/V sin err' respectively in these two cases, and we 
retrieve the integral equations for the even and odd prolate spheroidal 
functions of zero order. Note that when N = ±4, (25) reduces to the 
prolate spheroidal equation. 



VIII. NUMERICAL RESULTS 

A program for the IBM 7090 has been written to compute generalized 
prolate spheroidal functions using formulae (40) and (44). Trial values 
for the xjv.n were obtained from (34) and (55) and the recurrences 
(36)- (37) and (59). The method of Bouwkamp 6 was then used to cor- 
rect these estimates and obtain the d"'". Values of y N ,„ were obtained 
from (46) and these were converted to values of X by Xw.n = cyw.n • 

Fig. 1 shows plots of xw.n versus c. Fig. 2 gives the behavior of the 
first few Xat,„ . By definition of the labels, Xf.n+i ^ XN.n for N,n = 0, 1, 
. . . and if c > the inequality is strict. From Sturmian theory, it follows 
that Xff+i.n > XN, n • For the X's, one can show correspondingly that 
Xjv,„+i < Xiv,„ and \ N+ i, n < X.v.„ for N,n = 0, 1, . . . . The problem of 
ordering the X's and x's for all N and n appears to be a difficult one. Some 
values are listed on Table I. 

Figs. 3 and 4 show plots of <p N , n (x) versus x for N = 0,2, n = 0,1,2,3 
and c = 2,10. Values of the <ps,„ for a larger set of parameter values are 
given in Table II. Normalization is as in (86). 
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Fig. 2 — Curves of X N ,„ of (13) and (15) vs c. 
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Fig. 3 — Some generalized prolate spheroidal functions, <PN, n (x). 
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Fig. 4 — More generalized prolate spheroidal functions, *>Ar,„(j). 
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e 


Table I 

X 


— Numerical 1 

X 


Values of x-v,„ and 

c x 


X 






N = n = 






iV=l n 


= 




0.1 


7.5499895 - 1 


2.4968775 - 


3 


0.5 


3.9163765 4- 


9.4982658 - 


4 


0.5 


8.7434899 - 1 


6.0585348 - 


2 


1.0 


4.4119661 4- 


1.3986168 - 


2 


1.0 


1.2395933 + 


2.2111487 - 


1 


2.0 


6.3394615 4- 


1.6123183 - 


1 


1.5 


1.8225178 + 


4.2951906 - 


1 


3.0 


9.3427678 + 


4.8326866 - 


1 


2.0 


2.5857968 + 


6.2963045 - 


1 


4.0 


1.3086855 + 1 


7.8473505 - 


1 


3.0 


4.4622709 + 


8.8705036 - 


1 


5.0 


1.7170130 + 1 


9.3671678 - 


1 


4.0 


6.5208586 + 


9.7495117 - 


1 


6.0 


2.1310500 4- 1 


9.8534266 - 


1 


5.0 


8.5869176 + 


9.9534230 - 


1 


10.0 


3.7555900 4- 1 


9.9998314 - 


1 


10.0 


1.8690110 + 1 


9.9999957 - 


1 












N = n = 1 






N = 1 n 


= 1 




1 


9.2562398 4- 


1.0829815 - 


4 


1 


1.6255011 4- 1 


7.4672551 - 


7 


2 


1.0847476 + 1 


6.7214485 - 


3 


2 


1.7912353 + 1 


1.8549511 - 


4 


3 


1.3698728 + 1 


6.6745424 - 


2 


4 


2.4832293 4- 1 


3.8313651 - 


2 


4 


1.7898720 + 1 


2.6742780 - 


1 


5 


3.0401459 4- 1 


1.6818804 - 


1 


5 


2.3241561 + 1 


5.7877057 - 


1 


6 


3.7326440 4- 1 


4.2912557 - 


1 


(i 


2.9277622 + 1 


8.3060712 - 


1 


7 


4.5219234 4- 1 


7.1473948 - 


1 


7 


3.5550580 4 1 


9.4973850 - 


1 


8 


5.3565692 4- 1 


8.9618892 - 


1 


8 


4.1805821 + 1 


9.8782700 - 


1 


9 


6.1976089 + 1 


9.7041388 - 


1 


9 


4.7985976 + 1 


9.9738554 - 


1 


10 


7.0297509 + 1 


9.9279210 - 


1 


10 


5.4108072 4- 1 


9.9947801 - 


1 












N = n = 2 






JV=1 n 


= 2 




1 


2.5751488 + 1 


1.8834675 - 


9 


1 


3.6265101 + 1 


4.6976877 - 


12 


2 


2.6773866 + 1 


1.9235204 - 


6 


2 


3.7820310 4- 1 


1.9082396 - 


8 


5 


3.8241737 + 1 


1.6017987 - 


2 


5 


4.9160037 4- 1 


1.0305314 - 


3 


6 


4.4846367 4 1 


8.1254764 - 


2 


6 


5.5464880 4- 1 


8.3652966 - 


3 


7 


5.3021146 4- 1 


2.5847455 - 


1 


7 


6.3286568 4- 1 


4.4641026 - 


2 


8 


6.2527715 4- 1 


5.3544699 - 


1 


8 


7.2759605 + 1 


1.6080875 - 




9 


7.2854528 + 1 


7.8635574 - 


1 


9 


8.3789365 + 1 


3.9082845 - 




10 


8.3461406 + 1 


9.2600949 - 


1 


10 


9.5955815 4- 1 


6.6399691 - 




11 


9.4019226 + 1 


9.7915064 - 


1 


11 


1.0867089 4- 2 


8.6120123 - 




12 


1.0443896 + 2 


9.9484586 - 


1 


12 


1.2145589 4- 2 


9.5495434 - 




13 


1.1474313 4- 2 


9.9882732 - 


1 


13 


1.3409696 4- 2 


9.8759852 - 




14 


1.2496987 4- 2 


9.9974793 - 


1 


14 


1.4657506 + 2 


9.9692509 - 




15 


1.3514611 + 2 


9.9994820 - 


1 










Hi 


1.4528810 + 2 


9.9998984 - 


1 












N = n = 3 






N = 1 n 


= 3 




1 


4.9250694 + 1 


6.0066949 - 


15 


1 


6.4258409 4- 1 


7.6540787 - 


18 


2 


5.0761114 + 1 


9.8333952 - 


11 


2 


6.5789319 + 1 


4.9988893 - 


13 


5 


6.1688709 4 1 


3.5422330 - 


5 


5 


7.6749767 + 1 


1.1190662 - 


6 


7 


7.4995083 4- 1 


3.5278392 - 


3 


9 


1.0834214 + 2 


1.0298512 - 


2 


8 


8.3823340 4- 1 


2.0130790 - 


2 


10 


1.2001776 + 2 


4.5573797 - 


2 


9 


9.4336396 + 1 


8.2918248 - 


2 


11 


1.3350648 4- 2 


1.5017502 - 


1 


10 


1.0659367 + 2 


2.4212641 - 




12 


1.4872078 + 2 


3.5789707 - 


1 


11 


1.2034708 + 2 


4.9658387 - 




13 


1.6522672 4- 2 


6.1969392 - 


1 


12 


1.3504432 4- 2 


7.4660703 - 




14 


1.8237982 4- 2 


8.2818634 - 


1 


13 


1.5007176 + 2 


9.0244395 - 




15 


1.9961900 4- 2 


9.3866345 - 


1 


14 


1.6502439 4- 2 


9.6944048 - 




16 


2.1666412 + 2 


9.8150785 - 


1 


15 


1.7977291 + 2 


9.9165294 - 




17 


2.3347382 4- 2 


9.9500835 - 


1 


1(> 


1.9433894 4- 2 


9.9791376 - 
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Table I — 


- Continued 












N 


— 2 


n = 






.V = 


2 


n 


= 2 




1 


9.4970317 


+ o 


3.9517707 - 


4 


1 


4.9292042 


+ 


1 


1.0614368 - 


14 


2 


1.171091(5 


+ 1 


1.9088335 - 


2 


2 


5.0922802 


+ 


1 


1.7074109 - 


in 


3 


1.5291960 


+ 1 


1.3627864 - 


1 


5 


6.2566028 


+ 


1 


5.5396960 - 


5 


-1 


2.0048498 


+ 1 


4.0298411 - 


1 


7 


7.6509011 


+ 


1 


5.0165514 - 


3 


5 


2.5667098 


+ 1 


7.0317221 - 


1 


8 


8.5676381 


+ 


1 


2.7077307 - 


2 


fi 


3.1747966 


+ 1 


8.9417071 - 


1 


9 


9.6519700 


+ 


1 


1.0455606 - 


1 


7 


3.7952889 


+ 1 


9.7092917 - 


1 


10 


1 .0904728 


+ 


2 


2.8443961 - 


1 


8 


4.4125829 


+ 1 


9.9324606 - 


1 


11 


1.2295597 


+ 


2 


5.4626700 - 


1 


11) 


5.6324064 


+ 1 


9.9972026 - 


1 


12 


1.3768257 


+ 


2 


7.8216583 - 


1 












13 


1.5265435 


+ 


•> 


9.1919613 - 


1 












14 


1.6752489 


+ 


■> 


9.7528766 - 


1 












15 


1.8220066 


+ 


■> 


9.9334500 - 


1 




A' 


= 2 


n = 1 






iV = 


2 


n 


= 3 




1 


2.5333581 


+ 1 


4.1659113 - 


9 


1 


8.1275289 


+ 


1 


9.2115325 - 


21 


2 


2.7088321 


+ 1 


4.0530517 - 


6 


2 


8.2854667 


+ 


1 


2.3940665 - 


15 


fi 


3.9788041 


+ 1 


2.5787154 - 


2 


5 


9.4065006 


+ 


1 


3.2594658 - 


S 


fi 


4.6842565 


+ 1 


1.1592915 - 


1 


10 


1.3679420 


+ 


2 


5.8064778 - 


3 


7 


5.5371030 


+ 1 


3.2573165 - 


1 


11 


1.4973956 


+ 


2 


2.6682390 - 


2 


s 


6.5067655 


+ 1 


0.0754452 - 


1 


12 


1.6450899 


+ 


2 


9.5057491 - 


2 


9 


7.5425367 


+ 1 


8.3161633 - 


1 


13 


1.8113152 


+ 


2 


2.5380401 - 


1 


10 


8.5969115 


+ 1 


9.4455397 - 


1 


14 


1.9931616 


+ 


2 


4.9785798 - 


1 


ll 


9.6442775 


+ 1 


9.8484078 - 


1 


15 


2.1846670 


+ 


2 


7.3845601 - 


1 












16 


2.3793070 


+ 2 


8.9389256 - 


1 












17 


2.5727411 


+ 


2 


9.6462917 - 


1 












18 


2.7635393 


+ 


2 


9.8968676 - 


1 



Table II — Values of <p N ,n(x) 


,Y = » = 


X 


t = 1 r = 2 c = 5 


c = 10 



0.1 


4.74638 


— 


1 


5.55421 


— 


1 


9.15062 


— 


1 


1.31455 + 


0.2 


6.68776 


_ 


1 


7.74706 


— 


1 


1.22032 


+ 





1.62247 4- 


0.3 


8.14070 


— 


1 


9.27005 


— 


1 


1.35165 


+ 





1.57689 4- 


0.4 


9.31948 


— 


1 


1.03607 


+ 





1.35103 


+ 





1.30428 + 


0.5 


1.03044 


+ 





1.11011 


+ 





1.24626 


+ 


I) 


9.31637 - 1 


0.6 


1.11351 


+ 





1.15353 


+ 





1.06660 


+ 





5.70325 - 1 


0.7 


1.18341 


+ 





1 . 16921 


+ 





8.43474 


— 


1 


2.91331 - 1 


0.8 


1.24157 


+ 





1 . 15957 


+ 





0.07845 


— 


1 


1.17077 - 1 


0.9 


1.28896 


+ 





1 . 12695 


+ 





3.86969 


— 


1 


3.19741 - 2 


1.0 


1.32627 


+ 





1.07383 


+ 





2.01532 


— 


1 


3.00159 - 3 


1.1 


1.35405 


+ 





1.00285 


+ 





6.38588 


— 


2 


-1.09501 - 3 


1.2 


1.37278 


- 





9.1(1840 


— 




-2.24980 


— 


2 


4.23236 - 4 


1.3 


1.38285 


+ 





8.18791 


— 




-6.18395 


— 


2 


6.58696 - 4 


1.4 


1.38464 


+ 





7.11797 


— 




-(1.43066 


— 


2 


-2.21883 - 4 


1.5 


1.37853 


+ 





5.98995 


— 




-4.32142 


— 


2 


-5.95391 - 4 


1.6 


1.36489 


+ 





4.83499 


— 




-1.22744 


— 


2 


-1.21194 - 4 


1.7 


1.34410 


+ 





3.68328 


— 




1.68342 


— 


2 


4.20886 - 4 


1.8 


1.31655 


+ 





2.56332 


— 




3.61314 


— 


2 


3.79898 - 4 


1.9 


1.28264 


+ 





1.50130 


— 




4.20554 


— 


2 


-8.00341 - 5 
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Table II — Continued 



N = n = 



2.0 


1.24281 + 


5.20429 - 2 3.52268 - 


- 2 


-3.67577 - 


- 4 


2.1 


1.19748 + 


-3.59472 - 2 1.94309 - 


- 2 


-2.42458 - 


- 4 


2.2 


1.14713 + 


-1.12245 - 


1.30129 - 


- 4 


1.44876 - 


- 4 


2.3 


1.09221 + 


-1.75666 - 


-1.71274 - 


- 2 


3.08117 - 


- 4 


2.4 


1.03322 + 


-2.25459 - 


-2.80628 - 


- 2 


1.76068 - 


- 4 


2.5 


9.70655 - 1 


-2.61300 - 


-3.05323 - 


- 2 


-1.11428 - 


- 4 


2.6 


9.05010 - 1 


-2.83289 - 


-2.47697 - 


- 2 


-3.05239 - 


- 4 


2.7 


8.36800 - 1 


-2.91927 - 


-1.30522 - 


- 2 


-6.14275 - 


- 6 


2.8 


7.66537 - 1 


-2.88082 - 


1.07933 - 


- 3 


1.53741 - 


- 4 


2.9 


6.94735 - 1 


-2.72956 - 


1.38275 - 


- 2 


6.59394 - 


- 5 


3.0 


6.21906 - 1 


-2.48005 - 


I 2.20781 - 


- 2 


1.15331 - 


- 4 



N = n = 1 



c = 5 



c = 10 



0.1 


7.57682 




1 


7.49125 - 1 


7.51850 


- 1 


1.11517 + 


0.2 


1.00189 


+ o 


9.78062 - 1 


8.84348 


- 1 


9.64135 - 1 


0.3 


1.08562 


+ 





1.03485 + 


7.51864 


- 1 


2.39045 - 1 


0.4 


1.02660 


+ o 


9.38690 - 1 


4.08540 


- 1 


-6.68891 - 1 


0.5 


8.24694 




1 


6.94825 - 1 


-7.50154 


- 2 


-1.35733 + 


0.6 


4.76162 


— 


1 


3.08610 - 1 


-6.10130 


- 1 


-1.58691 + 


0.7 


-2.29984 


_ 


2 


-2.11235 - 1 


-1.10302 


+ 


-1.36404 4- 


0.8 


-6.75867 


_ 


1 


-8.51754 - 1 


-1.47075 + 


-8.87374 - 1 


0.9 


-1.48414 


+ 





-1.59571 + 


-1.65550 


+ 


-4.07419 - 1 


1.0 


-2.44790 


+ 





-2.41297 + 


-1.63388 


+ o 


-9.25172 - 2 


1.1 


-3.56556 


+ 





-3.30607 + 


-1.41956 


+ o 


2.63841 - 2 


1.2 


-4.83382 


+ o 


-4.22107 4- 


-1.05901 


+ o 


2.11679 - 2 


l.:} 


-6.24769 


+ 





-5.13836 + 


-6.20954 


— 1 


-1.29202 - 2 


1.4 


-7.80051 


+ 





-6.02868 + 


-1.82563 


— 1 


-2.10428 - 2 


1.5 


-9.48401 


+ 





-6.86308 + 


1.85537 


— 1 


-2.91179 - 3 


1.6 


-1.12884 


+ 


1 


-7.61390 + 


4.31871 


— 1 


1.47146 - 2 


1.7 


-1.32026 


4- 


1 


-8.25574 4- 


5.32289 


— 1 


1.36618 - 2 


1.8 


-1.52139 


+ 


1 


-8.76628 + 


4.92313 


— 1 


-1.33175 - 3 


1.9 


-1.73086 


+ 


1 


-9.12711 + 


3.43349 


— 1 


-1.28433 - 2 


2.0 


-1.94720 + 


1 


-9.32424 + 


1.34214 


— 1 


-9.24495 - 3 


2.1 


-2.16882 


+ 


1 


-9.34869 + 


-8.02848 


— 2 


2.21756 - 3 


2.2 


-2.39406 


+ 


1 


-9.19674 4- 


-2.51019 


— 1 


1.09934 - 2 


2.3 


-2.62120 + 


1 


-8.87007 + 


-3.44099 


— 1 


7.19680 - 3 


2.4 


-2.84846 


+ 


1 


-8.37588 + 


-3.46476 


— 1 


-2.21224 - 3 


2.5 


-3.07405 + 


1 


-7.72616 + 


-2.66729 


— 1 


-6.65003 - 3 


2.6 


-3.29612 


+ 


1 


-6.93789 + 


-1.31077 


— 1 


-7.12403 - 3 


2.7 


-3.51286 


+ 


1 


-6.03209 + 


2.39279 


— 2 


8.00041 - 4 


2.8 


-3.72245 


+ 


1 


-5.03334 + 


1.60443 


— 1 


8.27644 - 3 


2.9 


-3.92310 


+ 


1 


-3.96899 + 


2.48208 


— 1 


1.58822 - 3 


3.0 


-4.11308 + 


1 


-2.86785 + 


2.70290 


— 1 


-1.20582 - 3 







N = n = 2 










.T 


c = 1 


c = 2 


c = 5 




c = 10 




0.1 
0.2 


9.39351 - 1 
1.08208 + 


9.35161 - 1 
1.06262 4- 


8.84844 - 
9.01678 - 


- 1 

- 1 


9.05937 - 
5.16232 - 


1 

1 
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Table I] 

^Y = 


Q 


- Continued 










X 


c = 1 






c = 2 






c = 5 






c = 10 




0.3 


8.60349 


_ 


1 


8.22287 





1 


5.05390 


_ 


1 


-3.42126 


_ i 


0.4 


3.(14401 


— 


l 


2.90052 


— 


l 


-1.35880 


— 


1 


-9.82679 


— 1 


0.5 


-3.055(H) 


_ 


1 


-3.80932 


— 


l 


-7.81957 


— 


l 


-8.93819 


— 1 


0.0 


-9.56032 


— 


1 


-1.00914 


+ 


o 


- 1 . 10120 


+ 





-7.44970 


— 2 


0.7 


-1.32240 + 





-1.31928 


4- 


(i 


-1.03507 


4- 


i) 


9.85253 


— 1 


0.8 


-1.05832 


+ 





-9.83017 


— 


l 


-2.04525 


— 


1 


1.05901 


+ o 


0.1) 


2.05717 




1 


3.74908 


— 


l 


1.15147 


4- 





1.00911 


+ o 


1.0 


3.10217 


+ 





3.10040 


+ 


(i 


3.05038 


4 





9.70106 


— 1 


1.1 


8.22453 


4- 


1) 


7.79020 


+ 


(i 


5.10083 


+ 





2.03097 


— 1 


1.2 


1.01238 


4- 


1 


1.40703 


4- 


l 


7.09000 


4- 





-2.55403 


— 1 


1.3 


2.70035 -\ 


1 


2.41735 


+ 


l 


8.48753 


4- 





-2.08864 


— 1 


1.4 


4.34731 


+ 


1 


3.00100 


+ 


l 


9.03589 


+ 





-2.22585 


— 2 


1.5 


0.40001 


4- 


1 


5.22380 4 


l 


8.58154 


+ 





1.83653 


— 1 


1.6 


9.19010 


+ 


1 


7.11812 


+ 


l 


7.14579 


+ 





1.76198 


— 1 


1.7 


1.26331 


+ 


2 


9.34716 4 


l 


4.93120 + 





9 . 12279 


— 3 


l.S 


1.08874 


+ 


2 


1.19010 


4- 


2 


2.28768 


+ 





-1.36988 


— 1 


1.9 


2.20527 


+ 


2 


1.47564 


+ 


2 


-3.55354 


— 


1 


-1.39402 


— 1 


2.0 


2.82291 


4- 


2 


1.78757 


+ 


2 


-2.57167 


4- 


(J 


-1.18170 


— 2 


2.1 


3.55155 


+ 


2 


2.12083 


+ 


2 


-4.02160 


4- 





1 .05884 


— 1 


2.2 


4.40084 


4- 


2 


2.40904 


+ 


2 


-4.51700 


+ 





1 . 18832 


— 1 


2.3 


5.38002 


+ 


2 


2.82468 


+ 


2 


-4.05482 


+ 





1.72440 


_ 2 


2.1 


0.49779 


+ 


2 


3.17933 


4- 


2 


-2.81239 


+ 





-8.73616 


- 2 


2.5 


7.70218 


4- 


2 


3.52375 


+ 


2 


-1.10632 


+ 





-7.36407 


- 2 


2.0 


9.18038 


4- 


2 


3.84832 


+ 


2 


6.70904 


— 


1 


-2.63630 


- 2 


2.7 


1.07586 


4- 


3 


4 . 14324 


4- 


2 


2.10111 


4- 





3.13871 


- 2 


2.8 


1.25021 


4- 


:i 


4.39S92 


+ 


2 


3.05058 


4- 





9.47047 


- 2 


2.9 


1.44146 4- 


3 


4.00041 


+ 


2 


3.21902 


+ 





7.91312 


- 3 


3.0 


1.04988 


+ 


3 


4.75731 

.v = 


4- 

(I 


2 


2.66304 

3 


4- 





-5.44100 


_ 2 


X 


c = 1 






c = 2 






c = 5 






c = 10 




0.1 


1.04335 


+ 





1.03904 


4- 





1.00390 


+ 





8.10002 


- 1 


0.2 


9.41800 




1 


9.22225 


— 


1 


7.82101 


— 


1 


2.80879 


- 1 


0.3 


2.91804 


— 


1 


2.54987 


— 


1 


1.09411 


— 


2 


-5.85742 


- 1 


0.4 


-5.60503 


_ 


1 


-6.04220 


— 


1 


-8.19610 


— 


1 


-9.01249 


- 1 


0.5 


-1.10122 


+ 





-1.17035 


+ 





-1.15387 


4- 





-2.96465 


- 1 


0.G 


-1.04537 


+ 





-1.00493 


+ 





-0.55699 


— 


1 


7.15663 


- 1 


0.7 


-0.94405 




2 


3.87818 


— 


3 


5.00570 


— 


1 


1 . 12898 


+ 


0.8 


1.23039 


+ 





1.27745 


+ 


o 


1.45420 


4- 





2.99477 


- 1 


0.9 


1.17030 


+ 





1.12044 


4- 





0.85065 


— 


1 


-1.44767 


+ 


1.0 


-3.74103 


+ 





-3.74119 


+ 





-3.72277 


4- 





-2.98825 


+ 


1.1 


-1.94508 


+ 


1 


-1.80263 


+ 


1 


-1.36118 


4- 


1 


-3.23654 


4- 


1.2 


-5.51898 


4- 


1 


-5.13331 


+ 


1 


-3.01767 


+ 


1 


-1.98417 


+ 


1.3 


-1.24305 


+ 


2 


- 1 . 12500 


4- 


2 


-5.34231 


4- 


1 


-6.65948 


- 2 


1.1 


-2.45530 


+ 


2 


-2.15895 


+ 


2 


-8.18229 


4- 


1 


1.27315 


4- 


1.5 


-4.43417 


4- 


2 


-3.78480 


4- 


2 


-1.12269 


+ 


2 


1.31011 


+ 


1.6 


-7.50023 


+ 


2 


-6.20420 


+ 


2 


-1.40431 


4- 


2 


3.14806 


- 1 


1.7 


-1.20572 


+ 


3 


-9.64812 


4- 


2 


-1.01419 


+ 


2 


-7.44769 


- 1 


1.8 


-1.80039 


+ 


3 


-1.43724 


+ 


3 


-1.70733 


+ 


2 


-1.00264 


4- 


1.9 


-2.77451 


+ 


3 


-2.00513 


4- 


3 


-1.65204 


+ 


2 


-4.10710 


- 1 


2.0 


-4.02027 


4- 


3 


-2.87083 


+ 


3 


-1.43828 


+ 


2 


4.24001 


- 1 


2.1 


-5.08257 


4- 


3 


-3.! 1(101 12 


+ 


3 


-1.08200 


+ 


2 


8.04224 


- 1 
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Table II — Continued 

N = « = 3 



2.2 


-7.85998 + 3 


-5.16340 + 3 


-6.24706 4- 1 


4.53989 - 1 


2.3 


-1.06656 + 4 


-6.68939 + 3 


-1.27832 + 1 


-2.40332 - 1 


2.4 


-1.42276 + 4 


-8.49881 + 3 


3.37407 + 1 


-6.90334 - 1 


2.5 


-1.86902 + 4 


-1.06060 4- 4 


7.02754 + 1 


-3.30234 - 1 


2.6 


-2.42135 + 4 


-1.30186 4- 4 


9.15596 4- 1 


1.36373 - 1 


2.7 


-3.09743 4- 4 


-1.57357 + 4 


9.49226 4- 1 


2.58923 - 1 


2.8 


-3.91655 + 4 


-1.87469 + 4 


8.08334 + 1 


4.73302 - 1 


2.9 


-4.89962 + 4 


-2.20322 + 4 


5.29160 + 1 


-7.99752 - 3 


3.0 


-6.06911 + 4 


-2.55591 + 4 


1.70562 4- 1 


-4.71948 - 1 



N = \ « = 



c = 10 



0.1 

0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 
1.1 
1.2 
1.3 
1.4 
1.5 
1.6 
1.7 
1.8 
1.9 
2.0 
2.1 
2.2 
2.3 
2.4 
2.5 
2.6 
2.7 
2.8 
2.9 
3.0 



+ 

+ 

+ 

+ 

+ 



6.67799 - 2 
1.88413 
3.44707 
5.27637 
7.31901 
9.53337 
1.18837 
1.43380 
1.68661 
1.94398 
2.20321 
2.46172 + 
2.71702 + 
2.96675 + 
3.20863 + 
3.44053 + 
3.66041 + 
3.86640 4- 
4.05675 + 
4.22984 + 
4.38426 + 
4.51871 + 
4.63210 4- 
4.72350 4- 
4.79216 + 
4.83753 4- 
4.85923 4- 
4.85708 + 
4.83107 + 
4.78140 + 



7.82376 
2.19147 
3.96102 
5.96043 
8.08692 
1.02496 
1.23655 4- 
1.43586 4- 
1.61605 4- 
1.77112 4- 



1.89601) 
1.98689 
2.04079 
2.05612 
2.03244 
1.97051 
1.87221 
1.74049 
1.57927 
1.39322 
1.18771 
9.68543 
7.41814 
5.13706 
2.90272 
7.72847 
-1.19954 
-2.96762 
-4.49226 
-5.74265 



- 2 

- 1 

- 1 

- 1 

- 1 
4- 












+ 



+ 

+ 

4- 

4- 

+ 

4- 

+ 

4- 

4- 



1.75066 - 1 
4.70668 - 1 
7.93680 - 1 
1.08116 4- 
1.28498 4- 
1.37488 4- 
1.34097 + 
1.19359 4- 
9.60016 
6.78651 
3.91685 
1.37762 

-5.42516 

-1.69153 

-2.06531 

-1.79111 

-1.08794 

-2.14572 
5.85005 
1.12680 
1.31527 
1 . 15060 
7.15483 
1.46127 

-4.05001 

-8.05860 

-9.73081 

-8.86612 

-5.89386 

-1.69513 



3.92683 - 1 

9.77051 - 1 

1.44397 + 

1.62479 4- 

1.49159 4- 

1.13769 4- 

7.13592 

3.52409 

1.21545 

1 . 73276 
-6.46094 
-6.17021 

4.36117 

1.77204 
-2.50040 
-2.92887 

5.03281 

2.47164 

1.94403 
-5.40072 
-2.13282 
-1.41516 

6.26047 

1.92158 

7.95172 
-6.86022 
-7.63503 
-9.86627 

1.47764 

1.50302 









N = 1 n = 








X 


c = 1 




c = 2 


c = 5 




c = 10 


0.1 
0.2 


1.78797 - 
4.81623 - 


- 1 

- 1 


1.86209 - 1 
4.98465 - 1 


2.26456 - 
5.77313 - 


- 1 

- 1 


4.69291 - 1 
1.01396 4- 
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Table II — Continued 











N = 


l n = 


l 








X 


c -1 






c = 2 




c = 5 






f = 10 


0.3 


8.11322 




1 


8.30025 


- 1 


8.81124 


_ 


1 


1.10220 + 


0.4 


1.0!) 107 


+ 





1.09943 


+ o 


1.01214 


+ o 


5.09582 - 1 


0.5 


1.24408 


+ 


1) 


1.22371 


+ o 


8.87050 


— 




-3.54754 - 1 


0.6 


1 . 19004 


+ 





1 . 12484 


+ o 


4.82831 


— 




-1.21477 + 


0.7 


8.42470 




1 


7.29000 


- 1 


-1.03019 


— 




-1.61645 + 


0.8 


1.15744 


_ 


1 


-2.01055 


_ 2 


-9.50512 


— 




-1.44879 4- 


0.0 


-1.07772 


+ 





-1.19400 


+ o 


-1.70308 


+ 




-9.00520 - 1 


1.0 


-2.82510 


+ 





-2.81285 


+ o 


-2.43994 


+ 




-3.20420 - 1 


1.1 


-5.21224 


f 





-4.90037 


+ o 


-2.80233 


+ 




3.12085 - 2 


1.2 


-8.32175 


+ 





-7.45742 


+ o 


-2.94990 


+ 




1.08983 - 1 


1.3 


-1.22323 


+ 


1 


-1.04040 


+ 1 


-2.08284 


+ 




2.38728 - 2 


1.4 


-1.70172 


+ 


1 


-1.38798 


+ 1 


-2.10030 


+ 




-0.19490 - 2 


1.5 


-2.27431 


+ 


1 


-1.70444 


+ 1 


-1.32097 


+ 




-0.25904 - 2 


1.6 


-2.94003 


+ 


1 


-2.10792 


+ 1 


-4.02707 


— 




-1.89305 - 4 


1.7 


-3.72401 


1 


1 


-2.58893 


+ 1 


3.24280 


— 




5.11305 - 2 


l.S 


-4.01145 


+ 


1 


-3.01663 


+ 1 


9 . 15552 


— 




4.39191 - 2 


1.9 


-5.01051 


+ 


1 


-3.43921 


+ 1 


1.22972 


+ 




-3.52397 - 3 


2.0 


-0.72374 


+ 


1 


-3.84424 


+ 1 


1.24232 


+ 




-3.87077 - 2 


2.1 


-7.95192 


+ 


1 


-4.21912 


+ 1 


9.87598 


— 




-3.49013 - 2 


2.2 


-9.29459 


+ 


1 


— 4.55151 


+ 1 


5.48435 


— 




2.81473 - 3 


2.3 


-1.07501 


+ 


2 


-4.82978 


+ 1 


3.07803 


— 




3.09193 - 2 


2.4 


-1.23154 


+ 


2 


-5.04353 


+ 1 


-4.30799 


— 




3.14921 - 2 


2.5 


-1.39803 


+ 


2 


-5.18371 


+ 1 


-7.57370 


— 




-0.17803 - 4 


2.0 


-1.57573 


+ 


2 


-5.24330 


+ 1 


-8.84103 


— 




-2.95208 - 2 


2.7 


-1.70214 


+ 


■> 


-5.21739 


+ 1 


-8.00250 


— 




-1.28822 - 2 


2.8 


-1.95708 


+ 


2 


-5.10350 


+ 1 


-5.41499 


— 




-1.34150 - 3 


2.0 


-2.15901 


-1- 


2 


-4.90184 


+ 1 


-1.79910 


— 




0.07715 - 3 


3.0 


2.30S71 


+ 


2 


-4.01480 


+ 1 


1.95783 






2.03000 - 2 


* 


r = 1 






.V = 

r = 2 


= 1 ft = 


2 

c = 5 






c = 10 


0.1 


3.17572 




1 


3.23375 


- 1 


3.58452 


— 


1 


4.01327 - 1 


0.2 


7.89038 


— 


1 


7.98450 


- 1 


8.45034 


— 


1 


8.84405 - 1 


0.3 


1.13929 


+ 





1.13924 


+ o 


1.10027 


+ 





6.93291 - 1 


0.1 


1.10208 4 





1 . 13700 


+ o 


9.09831 


— 


1 


-8.13380 - 2 


0.5 


7.40307 


— 


1 


0.88194 


- 1 


2.60402 


— 


1 


-8.00192 - 1 


0.6 


-0.41935 


— 


2 


-1.40156 


- 1 


-0.18349 


— 


1 


-9.72981 - 1 


0.7 


-9.90S30 


— 


1 


-1.04501 


+ 


-1.27959 


+ 





-2.01002 - 1 


o.s 


-1.44413 


+ 





-1.42980 


+ o 


-1.15587 


+ 





1.02436 + 


0.9 


-4.43333 




1 


-3.00470 


- 1 


2.92752 


— 


l 


1.92471 + 


1.0 


3.40330 


+ 


3.46067 


+ 


3.41429 


+ 





1.94059 + 


1.1 


1.22701 


+ 


1 


1 . 10998 


+ 1 


8.22013 


+ 


(1 


1.12202 + 


1.2 


2.80105 


+ 


1 


2.03729 


+ 1 


1.43419 


+ 


1 


5.54149 - 2 


1.3 


5.57952 


+ 


1 


4.98022 + 1 


2.09872 


+ 


1 


-5.92374 - 1 


1 .1 


0.78715 


+ 


1 


8.45483 


+ 1 


2.71147 


+ 


1 


-5.45072 - 1 


1.5 


1.59705 


»- 


2 


1.33319 


+ 2 


3.15908 


+ 


1 


-5.23392 - 2 


Mi 


2. 10994 


+ 


2 


1.98850 


+ 2 


3.34313 


+ 


1 


3.78477 - 1 


1.7 


3.0030S 


+ 


2 


2.S3S10 


+ 2 


3.20100 


+ 


1 


4.00284 - 1 


l.S 


5.25100 


+ 


2 


3.90597 


+ 2 


2.72591 


! 


1 


7.08751 - 2 


1.9 


7.31713 


+ 


2 


5.21250 


+ 2 


1.90583 


-1 


1 


-2.08803 - 1 


2.0 


9.95303 


+ 


2 


0.77290 


+ 2 


1.02408 


+ 


1 


-3.08083 - 1 
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Table II — Continued 











N = 1 « = 


■■ 2 








X 


c = 1 






c = 2 


c = 5 






c = 10 


2.1 


1.32612 


+ 3 


8.59566 + 2 


4.17984 


_ 


1 


-9.67834 - 2 


2.2 


1.73484 


+ 


3 


1.06812 + 3 


-8.32684 


+ 





1.95130 - 1 


2.3 


2. 2331(> 


+ 


3 


1.30213 + 3 


-1.46613 


+ 


1 


2.53304 - 1 


2.4 


2.83330 


+ 


3 


1.55981 + 3 


-1.76896 


+ 


1 


1.16918 - 1 


2.5 


3.54832 


+ 3 


1.83831 + 3 


-1.71224 


+ 


1 


-9.81521 - 2 


2.0 


4.39138 + 3 


2.13379 + 3 


-1.33267 


+ 


1 


-2.49959 - 1 


2.7 


" 5.37625 + 


3 


2.44140 + 3 


-7.25150 


+ 





-5.82649 - 2 


2.8 


6.51088 


+ 


15 


2.75538 + 3 


-2.34402 


— 


1 


9.65074 - 2 


2.9 


7.82720 


+ 


3 


3.06921 + 3 


0.27542 


+ 





5.32562 - 2 


3.0 


9.32135 


+ 


3 


3.37554 + 3 
N = 1 n — 


1.10027 

: 3 


+ 


1 


1.36992 - 1 


X 


c = 1 






c = 2 


c = 5 






c = 10 


0.1 


4.70292 


_ 


1 


4.74792 - 1 


5.03663 


_ 


1 


5.60653 - 1 


0.2 


1.03040 


+ 





1.03944 + 


1.05094 


+ 





9.52509 - 1 


0.3 


1.15201 


+ 





1.13873 + 


1.02612 


+ 





4.75110 - 1 


0.4 


5.87098 


— 


1 


5.49119 - 1 


2.80849 


— 


1 


-5.29321 - 1 


0.5 


-4.34815 


— 


1 


-4.78943 - 1 


-7.50431 


— 


1 


-1.03082 + 


0.6 


-1.22557 


+ 





-1.23658 + 


-1.23879 


+ o 


-3.74929 - 1 


0.7 


-9.49593 


— 


1 


-9.03408 - 1 


-5.28507 


— 


1 


8.87232 - 1 


0.8 


5.95761 


— 


1 


6.56832 - 1 


1.04428 


+ 


1.26903 + 


0.9 


1.58393 


+ 





1.56393 + 


1.35005 


+ 





-3.80514 - 1 


1.0 


-3.99973 


+ 


I) 


-3.99876 + 


-3.98415 + 


-3.66768 + 


1.1 


-2.94218 


+ 


1 


-2.82967 + 1 


-2.13760 


+ 


1 


-6.63805 + 


1.2 


-9.93028 


+ 


1 


-9.31025 + 1 


-5.82859 


+ 


1 


-7.12958 + 


1.3 


-2.54877 


+ 


2 


-2.33198 + 2 


-1.21470 


+ 


2 


-4.48170 + 


1.1 


-5.60512 


+ 


2 


-5.00166 + 2 


-2.14807 


+ 


2 


-1.91044 - 1 


1.5 


-1.11157 


+ 


3 


-9.66261 + 2 


-3.37164 


+ 


2 


3.09018 + 


l.(i 


-2.04369 


+ 


3 


-1.72812 + 3 


-4.80966 


+ 


2 


3.47603 + 


1.7 


-3.54344 


+ 


3 


-2.91000 + 3 


-6.31861 


+ 


2 


1.24407 + 


1.8 


-5.86049 


+ 


3 


-4.66630 + 3 


-7.70049 


+ 


2 


-1.54042 + 


1.9 


-9.32112 


+ 3 


-7.18302 + 3 


-8.72930 


+ 


2 


-2.73125 + 


2.0 


-1.43431 


+ 


1 


-1.06778 + 4 


-9.19104 


+ 


2 


-1.43583 + 


2.1 


-2.14517 


+ 


I 


-1.53984 + 4 


-8.92696 


+ 


2 


7.02031 - 1 


2.2 


-3.12971 


+ 


1 


-2.16194 + 4 


-7.87269 


+ 


2 


2.14894 + 


2.3 


-1.46725 


+ 


4 


-2.96368 + 4 


-6.08384 


+ 


2 


1.44329 + 


2.4 


-6.25326 


+ 


4 


-3.97614 + 4 


-3.74068 + 


2 


-1.90745 - 1 


2.5 


-8.60136 


+ 


4 


-5.23076 + 4 


-1.12896 + 


2 


-1.17683 + 


2.0 


-1.16452 + 


5 


-6.75841 + 4 


1 .40286 


+ 


2 


-1.61055 + 


2.7 


-1.55406 


+ 


5 


-8.58790 + 4 


3.50030 


+ 


2 


-7.92157 - 2 


2.8 


-2.04672 


+ 


5 


-1.07450 + 5 


4.86619 


+ 


2 


1.29119 + 


2.9 


-2.66305 


+ 


5 


-1.32504 + 5 


5.32282 


+ 


2 


3.73378 - 1 


3.0 


-3.42633 


+ 


5 


-1.01182 + 5 

iV = 2 n = 


4.83375 


+ 2 


4.32692 - 1 


*• 


c = 1 






c = 2 


c = 5 






c = 10 


0.1 


8.11214 


_ 


3 


9.30928 - 3 


2.21477 





2 


8.00048 - 2 


0.2 


4.58035 


— 


2 


5.22724 - 2 


1.20067 


— 


1 


4.01321 - 1 
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Table II 

A' = 


: 2 


- Continued 

n = 










X 


c = 1 






c = 1 






c = 5 






c = 10 




0.3 


1.25827 





1 


1.42274 





1 


3.07992 


_ 


1 


9.01985 


- 1 


0.4 


2.57171 


— 


1 


2.87029 


— 


1 


5.71015 


— 


1 


1.38094 


+ o 


0.5 


4.46741 


_ 


1 


4.90296 


— 


1 


8.72575 


— 


1 


1.62971 


+ o 


0.6 


6.99877 


— 


1 


7.52394 


— 


1 


1.16320 -r 





1.54932 


+ o 


0.7 


1.02059 


+ 


(1 


1 .07049 


+ 





1.39112 


+ 





1 . 19392 


+ o 


O.S 


1.41171 


+ 





1 .43889 


+ 





1.51294 


+ 


() 


7.25292 


- 1 


0.9 


1.87493 


+ 





1.84932 


+ 





1.50259 


+ 





3 . 16554 


- 1 


1.0 


2.41089 


+ 





2.29129 


+ 





1.35656 


+ 





6.79148 


- 2 


1.1 


3.01924 


+ 





2.75254 


+ 


1) 


1.09462 


+ 





-2.02873 


_ 2 


1.2 


3.69863 


+ 


3.21944 


+ 


(J 


7.56039 


— 


1 


-1.46045 


_ 2 


1.3 


4.44678 


+ 


I) 


3.67758 


+ 





3.91913 


— 


1 


1 .04085 


_ o 


1.1 


5.26054 


+ 





4.11222 


+ 





5.55434 


— 


2 


1.52034 


- 2 


1.5 


6.13588 -4 


t) 


4.50887 


+ o 


-2.07467 


— 


1 


1.29918 


- 3 


1.6 


7.06802 


+ o 


4.85377 


+ o 


-3.67156 


— 


1 


-1.11917 


- 2 


1.7 


8.05139 


+ 





5.13443 


+ 





-4.13815 


— 


1 


-9.67213 


- 3 


I.S 


9.07977 


+ 





5.34008 


+ 





-3.58514 


— 


1 


1.55354 


- 3 


L.g 


1.01403 


+ 


1 


5.46205 


+ 





-2.29434 


— 


1 


9.62330 


- 3 


2.0 


1.12430 


+ 


1 


5.49406 


+ 





-6.51947 


— 


2 


6.50298 


- 3 


2.1 


1.23637 


+ 


1 


5.43250 


+ 





9.35338 


— 


2 


-2.03792 


- 3 


2.2 


1.34985 


+ 


1 


5.27657 


+ 





2.12220 


— 


\ 


-8.19124 


- 3 


2.3 


1.46391 


+ 


1 


5.02833 


+ 





2.68832 


— 


1 


-5.06774 


- 3 


2.4 


1.57769 


+ 


1 


4.69273 


+ 





2.57171 


— 


1 


1.94735 


- 3 


2.5 


1.69026 


+ 


1 


4.27721 


+ 





1.86697 


— 


1 


4.94597 


- 3 


2.G 


1.80072 


+ 


1 


3.79179 


+ o 


7.90086 


— 


2 


5.04030 


- 3 


2.7 


1.90816 


+ 


1 


3.24856 


+ 





-3.78705 


— 


2 


-7.03146 


- 4 


2.8 


2.01167 


+ 


1 


2.66139 


+ o 


-1.36199 


— 


T 


-6.15552 


- 3 


2.9 


2.11036 


+ 


1 


2.04550 


+ 





-1.94874 


— 


1 


-1.13042 


- 3 


3.0 


2.20337 


+ 


1 


1.41675 

N - 


+ o 

= 2 •* — 


-2.03218 
1 




1 


1.09442 


- 3 


X 


c = 1 






c = 2 






c = 5 






c = 10 




0.1 


3.01525 


_ 


2 


3.18266 





2 


4.36114 


_ 


2 


1.10168 


- 1 


0.2 


1.63412 


_ 


1 


1.71721 


— 


I 


2.27802 


— 


1 


5.09403 


- 1 


0.3 


4.17552 


_ 


1 


4.35481 


— 


1 


5.46215 


— 


1 


9.75159 


- 1 


0.4 


7.63468 


— 


1 


7.87547 


— 


1 


9.08386 


— 


1 


1 . 10090 


+ 


0.5 


1 . 12468 


+ 





1 . 14274 


+ 





1.16775 


+ 





6.37139 


- 1 


0.6 


1.37439 


+ 





1.36734 


+ 





1.16137 


+ 





-2.87398 


- 1 


0.7 


1.33292 


+ 


1) 


1.28261 


+ 


(1 


7.56169 


— 


1 


-1.22293 


+ o 


0.8 


7.66020 




1 


6.71043 


— 


1 


-1.09012 


— 


l 


-1.69212 


+ o 


0.9 


-6.15776 


_ 


1 


— 7.14755 


— 


1 


-1.39595 


+ 


i) 


-1.50113 


+ o 


1.0 


-3.15683 


+ 


(1 


-3.14077 


+ 


1) 


-2.90182 


+ 





-8.45332 


- 1 


1.1 


-7.25498 


+ 





-6.87847 


+ 


(1 


-4.57844 


+ 





-1.47679 


- 1 


1.2 


-1.33591 


+ 


1 


-1.21904 


+ 


1 


-5.97349 


+ 





2.35668 


- 1 


1.3 


-2.19659 


+ 


1 


-1.93150 


+ 


1 


-6.88402 


+ 





2.24248 


- 1 


1.1 


-3.36151 


+ 


1 


-2.84518 


+ 


1 


-7.11308 


+ 


3.10851 


- 3 


1.5 


-4.88847 


+ 


1 


-3.97466 


+ 


1 


-6.57336 


+ 





-1.65861 


- 1 


1.6 


-6.83847 


+ 


1 


-5.32791 


+ 


1 


-5.31243 


+ 





-1.46461 


- 1 


1.7 


-9.27502 


+ 


1 


-6.90514 


+ 


1 


-3.50808 


+ 





2.53438 


- 3 


1.8 


-1.22634 


+ 


2 


-8.69800 


+ 


1 


-1.43922 


+ 





1.23063 


- 1 


1.9 


-1.58698 


+ 


2 


-1.06891 


+ 


2 


5.66677 


— 


1 


1 . 16626 


- 1 


2.0 


-2.01606 


+ 


2 


-1.28518 


- 


2 


2.19548 


+ 





3.46739 


- 3 


2.1 


-2.52012 


+ 


2 


-1.51503 


+ 


2 


3.20680 


+ 





-9.50870 


- 2 
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Table II — Continued 











N = 


= 2 


n = 


1 






X 


c = 1 






c = 2 






c = 5 




c = 10 


2.2 


-3.10555 + 


2 


-1.75404 


+ 


2 


3.47937 4- 





-1.00215 - 1 


2.3 


-3.77844 


+ 


2 


-1.99703 


+ 


2 


3.03214 + 





-1.01271 - 2 


2.4 


-4.54453 


+ 


2 


-2.23824 


+ 


2 


2.01098 + 





7.86142 - 2 


2.5 


-5.40910 4 


2 


-2.47130 


+ 2 


0.84580 - 


1 


0.25020 - 2 


2.0 


-0.37087 


+ 


2 


-2.08987 


+ 


2 


-0.09192 - 


1 


1.80542 - 2 


2.7 


-7.45191 


+ 


2 


-2.88715 


+ 


2 


-1.70387 + 


-2.83514 - 2 


2.8 


-8.03750 + 


2 


-3.05073 


+ 


2 


-2.39029 + 





-8.09042 - 2 


2.9 


-9.93033 


+ 


2 


-3.19203 


+ 


2 


-2.45263 + 





-0.05093 - 3 


3.0 


-1.13498 


+ 


3 


-3.28923 

N = 


+ 

s 2 


2 


-1.97549 + 

2 





4.94337 - 2 


X 


c = 1 






c = 2 






c = 5 




f = 10 


0.1 


0.92988 





2 


7.12073 


_ 


2 


8.51196 - 


2 


1.31002 - 1 


0.2 


3.52490 


— 


1 


3.00980 


— 


1 


4.18221 - 


1 


5.07004 - 1 


0.3 


8.01544 


— 


1 


8.14059 


— 


1 


8.93117 - 


1 


9.45350 - 1 


0.4 


1.20740 


+ 





1.21237 


+ 





1.21221 + 





7.07342 - 1 


0.5 


1.27474 


+ 





1.25311 


+ 





1.05255 4- 





-3.70902 - 2 


0.0 


7.57995 


— 


1 


7.02242 


— 


1 


2.84085 - 


1 


-9.13190 - 1 


0.7 


-3.27092 


— 


1 


-3.94312 


— 


1 


-8.22800 - 


1 


-1.05173 4- 


0.8 


-1.41304 


+ 


-1.43709 


+ 





-1.49250 + 





-7.12440 - 2 


0.9 


-1.01239 


+ 





-9.55450 


— 


1 


-4.82515 - 


1 


1.57998 4- 


1.0 


3.74003 


+ 





3.73507 


+ 


3.09332 + 





2.80770 4- 


1.1 


1.70108 + 


1 


1.08770 


+ 


1 


1.23963 4- 


1 


2.91420 4- 


1.2 


4.78738 + 


1 


4 . 45904 


+ 


1 


2.64247 + 


1 


1.00265 4- 


1.3 


1.04979 


+ 


2 


9.51393 


+ 


1 


4.50130 + 


1 


-0.78930 - 2 


1.4 


2.03280 


+ 


2 


1.79119 


+ 


2 


0.85921 + 


1 


-1.18788 4- 


1.5 


3.01809 


+ 


2 


3.09537 


+ 


2 


9.27854 + 


1 


-1.12986 4- 


1.6 


0.05080 


+ 


2 


5.01708 


+ 


2 


1.14715 + 


2 


-2.05881 - 1 


1.7 


9.03954 


+ 


2 


7.73357 


+ 


2 


1.30538 + 


2 


7.08003 - 1 


1.8 


1.47041 


+ 


3 


1.14308 + 


3 


1.30797 + 


2 


8.79453 - 1 


l.i) 


2.18840 


+ 


3 


1.03343 


+ 


3 


1.31147 + 


2 


3.17527 - 1 


2.0 


3.15489 


+ 


3 


2.20396 4 


3 


1.12984 4- 


2 


-4.11104 - 1 


2.1 


4.44000 


+ 


3 


3.05047 


+ 


3 


8.37540 + 


1 


-7.13208 - 1 


2.2 


0.11858 


+ 


3 


4.03107 


+ 


3 


4.08923 4- 


1 


-3.71099 - 1 


2.3 


8.27598 


+ 


3 


5.20570 + 


3 


7.34290 4- 


2.40288 - 1 


2.4 


1 . 10092 


+ 


4 


0.59547 


+ 


3 


-2.92519 + 


1 


0.17233 - 1 


2.5 


1.44270 


+ 


4 


8.21064 


4- 


3 


-5.75607 4- 


1 


2.78319 - 1 


2.0 


1.80502 


+ 


4 


1.00505 + 


4 


-7.35509 4- 


1 


-1.40009 - 1 


2.7 


2.38122 


+ 


4 


1.21320 


+ 


4 


-7.52754 4- 


1 


-2.33031 - 1 


2.8 


3.00582 


+ 


1 


1.44287 


+ 


4 


-0.32751 4- 


1 


-4.05877 - 1 


2.9 


3.75458 + 


4 


1.09309 


+ 


4 


-4.05587 4- 


1 


1.20902 - 2 


3.0 


4.04442 + 


4 


1.90136 

A T = 


+ 

1 


4 


-1.18748 4- 
3 


1 


4.27288 - 1 


X 


r = l 






c = 2 






c = 5 




c = 10 


0.1 


1.20908 





1 


1.29109 





1 


1.44081 - 


1 


1.94100 - 1 


0.2 


5.91772 


— 


1 


5.99200 


— 


1 


0.48002 - 


1 


7.75163 - 1 


0.3 


1.13308 -f 





1 . 13787 


+ 





1.15991 + 





1.00828 4- 
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Table II — Continued 

S m 2 n = 3 






■ 


X 

0.4 


c = 1 

1 . 22242 


+ 


() 


r = 2 

1.20735 + 


c = 5 

1.08268 


+ 


c = 10 
4.42410 - 1 


0.5 


5.02125 




1 


4.61442 - 1 


1.72834 


— 


1 


-7.08002 - 1 


0.6 


-7.28650 


— 


1 


-7.67417 - 1 


-9.99672 


— 


1 


-1.12281 4- 


0.7 


-1.34994 


+ 





-1.33947 + 


-1.20442 


+ 


-6.35200 - 2 


0.8 


-1.55815 




1 


-9.78895 - 2 


3.13440 


— 


1 


1.38983 + 


0.0 


1.72855 


+ 


I) 


1.73119 + 


1.69507 


+ 


6.83573 - 1 


1.0 


-4.24197 


+ 





-4.23990 + 


-4.22216 


4- 





-4.05639 4- 


1.1 


-4.29825 


+ 


1 


-4.14973 + 1 


-3.22563 


-!- 


l 


-1.17483 4- 1 


1.2 


-1.69966 


+ 


2 


-1.60436 + 2 


-1.05765 


+ 


2 


-1.81766 4- 1 


L.3 


-4.92189 


+ 


2 


-4.54612 + 2 


-2.55016 


+ 


2 


-1.86221 4- 1 


1.1 


-1.19761 


+ 


:s 


-1.08177 4- 3 


-5.12503 


4- 


2 


-1.15388 4- 1 


1.5 


-2.59483 


+ 


3 


-2.28970 4- 3 


-9.05339 


+ 


2 


-3.41475 - 1 


1.6 


-5.16494 


+ 


3 


-4.44652 4- 3 


-1.44608 


4- 


3 


8.52760 + 


1.7 


-9.62735 


+ 


3 


-8.07479 4- 3 


-2.12354 


+ 


3 


1.01259 + 1 


1.8 


-1.70216 


+ 


1 


-1.38880 4- 4 


-2.89670 


+ 


3 


4.31901 + 


1.9 


-2.88067 


+ 


4 


-2.28281 4- 4 


-3.69313 


+ 


3 


-3.63918 + 


2.0 


-4.69805 


+ 


1 


-3.61016 4- 4 


-4.41450 


+ 


3 


-7.22580 4- 


2.1 


-7.42193 . 


+ 


4 


-5.52126 4- 4 


-4.94883 


+ 


3 


-5.31249 + 


2.2 


-1.14040 4 


5 


-8.19885 4- 4 


-5.18833 


+ 


3 


8.57079 - 1 


2.3 


-1.70988 


+ 


5 


-1.18598 + 5 


-5.05014 


+ 


3 


5.14940 + 


2.4 


-2.50841 


+ 


5 


-1.67557 + 5 


-4.49592 


4- 


3 


5.60742 + 


2.5 


-3.60853 


+ 


5 


-2.31714 + 5 


-3.54624 


+ 


3 


5.87060 - 1 


2.0 


-5.10005 


+ 


5 


-3.14232 + 5 


-2.28432 


+ 


3 


-4.18090 + 


2.7 


-7.09289 


+ 


5 


-4.18531 4- 5 


-8.50402 


+ 


2 


-2.39443 4- 


2.8 


-9.72020 


+ 


5 


-5.48229 + 5 


5.76599 


+ 


2 


-1.90159 4- 


2.9 


-1.31415 


+ 


c. 


-7.07072 + 5 


1.80814 


+ 


3 


6.87950 - 1 


3.0 


-1.75464 


+ 


6 


-8.98702 + 5 


2.67521 


+ 


3 


4.80391 4- 



I am indebted to Airs. E. Sonnenblick for programming and carrying 
out the computations reported here. 



APPENDIX A 



.4 Perturbation Scheme 

We treat briefly the following problem. Eigenfunctions u n and eigen- 
values X„ of an operator L are assumed known. That is, we have 



Lu„ + X„w„ = 0, a = 0, 1,2, ... . 



(114) 



It is desired to find eigenfunctions yp n and eigenvalues x». of the perturbed 
equation 



(L - eM)^ 4- x^ = 0. 



(115) 
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It is assumed that the u H satisfy the boundary condition to be imposed 
on the \p's and that the u n are complete in some appropriate sense. We 
proceed further in a purely formal maimer. 
Substitute the series 

<A„ = Un + 2 e'Qi (116) 

/=i 

Xn = A„ + E«'«i (117) 

;=1 

into (115). Here in the notation we have suppressed the dependence of 
the Qj and o> on n. By equating to zero the coefficients of distinct powers 
of e, we find 

Lw„ + X„m„ = (118) 

LQj + \,,Qj = MQ,_! - J2 a k Qj- k , (119) 

A = l 

3 = 1, 2, . . . 

where we define Q = u„ . 

Now it frequently happens that the operator M is such that Mm,, can 
be expressed as a finite linear combination of the w's with constant coef- 
ficients. We assume this to be the case and write 

i 

<=-< (120) 

n = 0,1,2, ... . 

Here a is a positive integer, / a nonnegative integer, and the superscript 
i on 7 is not an exponent, but a label. 

If the u n are linearly independent, formal solution of system (119) 
is now straightforward. Set 

n 

Qj = 2 A k 'u n+ka 

AJ=0 (121) 

3 = 1, 2, ... . 

The A's of course depend on n, but for simplicity we have suppressed 
this fact in the notation. Again the superscript is a label, not an ex- 
ponent. Substitute (121) and (120) into (119). Setting the coefficient 
of u n equal to zero in the resultant expression yields 
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aj= E^V-T,,-^', j= 1,2 (122) 

k — i 

Requiring the coefficient of u, l+ma to vanish gives 

i _ ' _ 

(X n+mtt - K)Aj = Y,a k AJ k - 2 -A—* 7«+(«-*)« ,,__. 

*— i fc — ' (123) 

m = -./'/, -./7+i, •■•,.//; 3 -1,2,.... 
Here we have adopted the conventions 

• U J - 
if either 

| /,• | > jZ, or a/.' < -n, or fc = and j = 1, 2, ... 
A„" = 1, .4," = 0, fc * 0, a D = 0. 

Equations (122) and (123) together with these conventions permit 
successive determination of the a's and A's. The case / = 1 occurs fre- 
quently. The first few coefficients for this case are given below where 
we have set 



A_x ] = -h- ay ,r 
Ai i = -w 

n ., = -[/j„7, 1 i 7«+« _i + h—y n~ l y »-«'] 

« 2 7 i -I -1 

.4_ 2 = «-2a/'-a7.« 7«-a 

•> . 2 —If I' i "i 

A_r = (h- a yyn [-y„ + Yn-«] 

i 2 / ; \2 It , l>| 

Ax = (ha)y, l \-y« + y„+a I 

A 2 2 = /l2a/i a 7« 1 7-,+« 1 

a ;i - (A„) 2 7» I 7» + «" 1 (-7« + 7« + i°) + (ft-a)V , 7«-« 1 (-7n + 7»-a°) 
A, 8 = h a [Ai(y» ~ 7.+-') + Mi 1 - t-W'Ai"] 
A_x 8 = fe_ a [A_ a 2 ( 7B - 7—°) + M-i 2 - yn-*a l Aj\ 

< :\ -lit 3 • 

«4 = Al7„ +a + -l-l7»-a ■ 

More generally for this case one finds 
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A±/ = (-1)' II A*****-*).* 1 , ./' =1,2,... 

il ±0 V = ^^j" £ h ±ka [ ai - y ll±ka \ j = 2, 3, . . . 

fc=l 

;'-2 

'^±o-2) = ^±(j-2) 2-, h±i<> (124) 

/=1 
X [^2 + ^±(J+l)a7n±(J+l)a 7n±/o 

+ (ai — y„±i a °) 2 A±*«(ai - 7n±*«°)L 



I 

j - 3, 4, . . . . 



APPENDIX B 



(125) 



Evaluation of an Integral 

We establish here the formula (43). Let 

•'o 

- f Ky(xy)T NtR {y)dy 

Jo 
on using the notation of (21). Then 

= f o T N , n {y)\xYK N "{xy) + 2xyK N '(xy) + (*' - x)/^^)]*j (126) 

= f T X M\(-W ~ I + W 2 + z 2 

- x)K N (xy) + 2xyK N '(xy)]dy 

by (23). Here primes denote differentiation with respect to the argu- 
ment indicated. 
Now 

dy (1 - y) —&y— + [TIT + V 7 '- W = ° 
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with x given by (20). Multiply this equation by K N {xy) and integrate 
from zero to one. There results 

■ t '-("'> [i (1 - " !) tf + (Mr 1 + *) r '-] * 
- Jf r "-<»>[J (1 - ^^ +(*^+ *)*<«»>]* 

where we have integrated by parts and made use of the fact that 
7vjv(0) = 7', v .„(0) = 0. Carrying out the indicated differentiation, we 
find 

= £ 7'.v.„(//) [(1 - ;yV/vV'(.ry) - 2xyK l /{xy) 

+ (i^i^ + x) Kjr(av)] 4 

= - f ' 7',.„(//)[(-.«-V " 1 + # 2 + * - x)/v.v(.r//) 
•'ii 

+ 2xyK ff ' (xy)]dy. 

Equations (120) and (127) give 



dy (127) 



[^»+>»E+{^-(^ +Sfc+ s)( y + Ste+ i)}! 



n,»(.r) = 0. 



The only solution of this equation that vanishes for x = is 

_ . . . Jn+1i+\ \x) 

F...W = k — ^=- . 
Using ( 125), wo then have 

fc JiNw(») = f j N ( xy )^/^Ts.n(y)dy. (128) 

To determine fe, we have only to compare the coefficients of x 
on both sides of (128). In this way we find 

k 



2»+*»+ i T{N + 2n + 2) (129) 

= ( ~ 1} " f y s+ - n+i T Nn (y)<ly. 

2 N +°"T(N + n + l)n! in J 
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The integral here can be evaluated by using (27), (28) and known 
properties of the Jacobi polynomials. We have 

f y N+2n+i TsAy)dy 

Jo 

= (" t N T I! y 2N+2n+1 Pn ilffi) (l - 2y 2 )dy 

= 2 _ (N+I1+2) /« + NX- 1 £ (1 _ w) , V(1 _ u) np n i».o) (u)du 

Now the coefficient of u" in (1 — w)"is ( — 1)" and the coefficient of 
u n in P n ( ' v ' 0, (u) is( 2n + N J I T [Ref. 8, Vol. II, p. 169, Eq. (5)], 



so 



(1 - u) " = / 4rrl\ p " ( ' v,0)( " ) + S a * p * ( * ,0)(m) ' 
v »> ) 



It follows, then, that 

f (1 - M)"(l - u) n P n {Nfi) ( U )du 

J (1 - u) P n (u)P„ (u)du 



2/i + N 
n 



(-l)»2»+*+ 1 



( 2n : ") 



(2n + iV + 1 



where we have used the orthogonality of the Jacobi polynomials and 
the known normalization integral [see Ref. 5, page 68, Eq. (4.3.3), for 
example]. Combining these results, (129) yields 

and together with (128) this establishes (43). 

REFERENCES 

1. Slepian, D., and Pollak, H. O., Prolate Spheroidal Wave Functions, Fourier 

Analysis and Uncertainty— I, B. S.T.J. , 40, Jan., 1901, pp. 43-64. 

2. Landau, H. J., and Pollak, H. O., Prolate Spheroidal Wave Functions, Fourier 

Analysis and Uncertainty — II, B.S.T.J., 40, January, 1961, pp. 65-84. 



PROLATE SPHEROIDAL WAVE FUNCTIONS — IV 3057 

3 Fox, A. G., and Li, Thieve, Resonant Modes in a Maser Interferometer, 

B.S.T.J., 40, March, UMii, pp. 453-488. 

4 Riesz, F., and Sz.-Nagy, B., Functional Analysis, Frederick Ingar Co., New- 

York, 1955. , n „ . U1 . 

5. Szegci, G., Orthogonal Polynomials, Am. Math. Soc. Colloquium Publication 

23, New York, 1959. 
(i Bouwkamp, C. J., On Spheroidal Wave Functions of Order Zero, J. Math. 

Phys.,26, 1957, pp. 79-92. 

7. Slepian, D., Some Asymptotic Expansions for Prolate Spheroidal Wave func- 

tions, to be published. 

8. Krdelyi, A., Higher Transcendental Functions, McGraw-Hill Book Co., New- 

York, 1953. . , , 

9. Kogelnik, H., Modes in Optical Resonators, Section 5.1, chapter in Advances 

in Lasers, ed. Levine, A. K., Dekker Publishers, New- York — to appear. 

10. Slepian, D., Analytic Solution of Two Apodization Problems, to be published. 

11. Heurtley, J. C, Hvperspheroidal Functions — Optical Resonators with Cir- 

cular Mirrors, to'appcar in Proceedings of the Polytechnic Inst, of Brooklyn 
Symposium on Quasi -Optics. June 8-10, 1904. 

12. Born, M., and Wolf, E., Principles of Optics, Pergamon Press, 1959. 



